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Ph D 
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ANALYSIS OE THE STEADY STATE AND DYNAMIC PERFORMANCE 
OE VOLTAGE CONTROLLED INDUCTION MOTOR DRIVES 

This thesis is devoted to the analysis and simulation of 
voltage controlled induction motor drives using thyristors 
Both open-loop and closed-loop (with speed signal feedback) 
operations are considered Analytical results are backed by 
experimental investigations 

It is a well established fact that the speed control of 
squirrel cage induction motors of small and medium sizes with 
fan-type loads is feasible, by voltage control, over reasonable 
ranges Although the voltage controlled induction motor (VCIM) 
drives are inefficient compared to inverter-fed and cycloconver- 
ter drives, the main reasons for its application are due to its 
reliability, simplicity and low cost Phase controlled thyri- 
stors can be effectively used for voltage control applications 
Although many schemes are possible, a wye-connected induction 
motor with a pair of back-to-back connected thyristors m each 
line is usually preferred for speed control applications# 

The use of VCIM drive for speed oontrol applications re- 
quires the analysis of the system response during steady state 


and transient conditions In this context the main contributions 
of the thesis can be summarized as follows 

1 ) The application of boundary value approach (BVA) to the 
steady state analysis using digital simulation 

2) The development of hybrid models of induction motors sui- 
table for transient and steady state analysis 

3) Simplified, yet accurate, evaluation of the steady state 
performance of VCIM from the solution of nonlinear algebraic 
equations 

4) Development of a linearized state space model of the VOIM 
and its application for predicting the stability of the 
VOIM drive from eigenvalue analysis 

5) Development of transfer functions and the design and testing 
of closed-loop control system for the speed control 

A brief description of the work reported in this thesis is 
given below 

The BVA, which has been used earlier for the analysis of 
invert er-fed induction motor drives, is extended to the steady 
state analysis and simulation of the following cases, 

1) wye-connected induction motor (isolated neutral) with 

a pair of back-to-back connected thyristors in eaoh line, 

2) induction motor with delta connected thyristors at the 
open star point, 

3) wye -connected motor with a pair of back-to-back connected 
thyristor-diode combination in each line. 



BVA makes use of boundary relationships and superposition 
principles for the solution The main features of this approach 
are 

a) the machine equations are invariant while the different 
modes of operation are accounted for by altering the 
forcing voltages, 

b) formulation becomes general and simple 

Analytical and experimental results and a comparative per 
formance of the drives with different thyristor configurations 
are presented 

Two hybrid models of an induction motor, that are well 
suited for the digital simulation of machines with terminal 
constraints, are developed 

The special features of the models are 

a) the identity of the motor terminals either on the stator 
side or on the rotor side is retained 

b) time dependence of the inductance matrix is eliminated 

c) the models are represented by simple equivalent circuits 
so that formulation of the state equations for any ter- 
minal constraint is possible with direct application of 
network theory 

The following examples are considered for the steady state 
and transient analysis using one of the hybrid models 



1) a wye- connected induction motor with a pair of back-to- 
back connected thyristors m each line, 

2) a delta connected induction motor with a pair of back-to- 
back connected thyristors in each line, 

A comparative study of the performance of these drives are 
present ed 

A simplified analysis of a VCIM drive, during steady state, 
is described in which , the motor performance is directly related 
to the control variables This relationship can be effectively 
used for the design of speed controllers during closed-loop 
operation Another feature of this simplified analysis is that 
the steady state operating variables (phase angle <fr, fundamen- 
tal components of stator currents and rotor flux linkages and 
average torque) can be predicted accurately by solving three 
nonlinear algebraic equations for any given value of triggering 
angle, slip and source voltage The assumption made m the 
simplified analysis is that the rotor flux linkages referred 
to a synchronously revolving reference frame are constants 
during constant speed operation Validity of this assumption is 
confirmed analytically and justified by simulation results 

Using a linearized model, the stability analysis of the 
VCIM drive is performed Prom the eigenvalue analysis, it is 
observed that certain category of machines may experience 
instability at near maximum speed 



Using the linear model of the VCIM drive and its associated 
transfer functions, a speed controller for closed -loop operation 
is designed Classical control theory is used for the design 
of the controller Dynamic response of the closed-loop system 
due to perturbations in the load torque and reference input is 
determined both analytically and experimentally 



CHAPTER 1 


INTRODUCTION 


1 1 GENERAL 

The speed control of induction motor drives has been 
revolutionized with the advent of thyristors in 1957 Thy- 
ristors capable of switching several hundred amperes at quite 
high voltages provided a flexible and reliable method of 
power control. The impact of these developments was first 
felt in the area of power conversion - repla cement of bulky 
and costly mercury are rectifiers, thyratrons and motor- 
generator sets by thyristors Until the middle of sixties* 
the d.c machine in spite of its disadvantages - costly and 
regular maintenance - occupied the dominant role in the 
field of variable spbedcdriiTas The robust, cheap, rugged 
and practically maintenance free squirrel cage induction 
motor was mostly restricted to applications involving con- 
stant speed operation, due to its inherent constant speed 
characteristics Now, with large scale application of thy- 
ristors, squirrel cage induction motor drives are increasingly 
replacing the conventional d.c motors m speed control appli- 
cations 

Many thyristor controlled drives have been developed 
in the past two decides However, exact analysis and design 
of many of these systems are still being developed due to the 
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complexity of the problem introduced by the switching elements 
(thyristors) The currents and voltages are non-smusoidal 
and the forcing voltages are difficult to determine m many 
instances In this connection, modern digital computers have 
come to the help of designers and researchers to a great 
extent Computer based analysis and design hare boen developed 
for many complex systems and this topic has been an active 
subject of research Although conventional models of the 
induction machine have been extensively used for the analysis, 
many attempts have been mado to develop new digital models of 
the machine to cope with the new situation created by the 
application of thyristors 

The two common techniques employed for the speed control 
of squirrel cage induction motors are 

a) stator voltage control 

b) stator voltage and frequency control 

1.2 STATOR VOLTAGE CONTROL 

Simple and most reliable method of speed control of 
squirrel cage induction motor with fan-type load is by the 
use of stator voltage control Since the torque developed 
m an induction motor is proportional to the square of the 
r.m.s voltage applied to the motor, the torque-speed chara- 
cteristics of the motor (Eig 1 1) can be adjusted by varying 
the input voltage Conventionally variable voltage is ob- 
tained by using vanacs and tap changing transformers The 




same objective is achieved by using phase controlled thy- 
ristors 11-3] in -which the triggering angle (a) is varied to 
obtain variable voltage The phase controlled thyristors 
are turned off at the natural current zero and hence no forced 
commutation is required Speed control of induction motors 
fed from a fixed frequency supply is inevitably an inefficient 
process The slippage power ((« -a> r )T) - difference between 
the stator input power (w T) and the rotor output power (« T) - 
is dissipated m the rotor as heat { 4] The motor may ex- 
perience excessive heating at lower speeds and special care 
is to be taken while selecting a motor for voltage controlled 
applications 

Although the voltage control scheme for the speed control 
of induction motors is an inefficient approach, this is em- 
ployed for many industrial applications of small and medium 
type induction motors with nonlinear load characteristics 
The reason for its popularity is its simplicity and reliability 
The cost of a voltage controlled drive is very much less com- 
pared to frequency controlled drives 

1 3 STATOR VOLTAGE AND FREQUENCY CONTROL 

The most efficient method of speed control of squirrel 
cage induction motors is by means of variable frequency con- 
trol A variable frequency control is normally accompanied by 
voltage control as well so that the air-gap flux remains con- 
stant throughout the operating range Thyristor inverters f 51 
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are now commonly used for the variable frequency, variable 
voltage control A variable d c voltage, obtained by means 
of a controlled rectifier, is fed to an inverter to produce 
variable frequency and variable magnitude a c voltages The 
static inverters are designed with a variety of circuit conf_- 
gurations [6,7 1 and employ forced commutation methods for 
turning off of the thyristors Basically there are two types 
of inverters - voltage source inverter and current source 
inverter [8,91 

Variable frequency can also be obtained using cyclo- 
converters [101 A cycloconverter, which converts alter- 
nating current at one frequency to alternating current at 
another frequency, employs line or load commutation for 
turning off of thyristors Since, there 1 ° no d c link m 
cycloconverters, these are more efficient than inverter-fed 
syst ems 

Although cycloconverter and inverter-fed drives are 
more efficient, they are very expensive and hence applicable 
mainly for large machines 

1 4 REVIEW 

Before discussing the work reported m the thesis a 
brief literature survey, pertaining to the analysis of vol- 
tage controlled induction motor drive, is undertaken 
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1 4 1 Steady State Analysis 

Shepherd and Stanway [11] reported an approximate solu- 
tion for a delta connected machine with back -to -hack connected 
thyristors, m series with each phase winding An exact solu- 
tion was attempted by Shepherd [12] m which the motor was 
represented by five nonlinear differential equations in terms 
of winding currents and time dependent inductance coefficients 
Possible approaches to the solution are discussed but no solu- 
tion is attempted 

Takeuchi [131 proposed a $ -function method to analyse 
a wye-connect ed induction motor with three pairs of back-to- 
back thyristors, connected in delta, at the open neutral 
point The average torque, voltages and currents are appro- 
ximated as an infinite series of Pouner components and the 
computation of the resulting equations is tedious Since the 
back-emf across the open circuited phase is neglected, accu- 
racy of the results will be poor for large delay angles, when 
the motor is near maximum speed 

lipo [14] presented a steady state analysis, using state 
variable technique, of a wye-connected induction motor with 
thyristor voltage control Closed form expressions for the 
steady state initial vector are derived which lead to the 
steady state solution The different modes of operation 
during voltage control are taken into 'account by changing the 
machine equations Analysis of a similar scheme has been 


presented by Ramamoorty and Samek [15J and results obtained 
by state variable method and harmonic equivalent circuit are 
compared The basic idea behind harmonic equivalent circuit 
method is to determine the Fourier components of the stator 
voltage and the resulting current components using steady 
state harmonic equivalent circuits The back-emf and phase 
angle <J> were obtained by an iterative procedure 

Bedford and Rene f 1 6 1 described a time domain analysis 
of a wye-connected motor using a-@-o transformation Ilango 
and Ramamoorty [17] presented the analysis of a wye-connected 
four wire induction motor The mathematical model used for 
the analysis retains the identity of the stator currents for 
motors, with neutral connected or unconnected 

Lipo's [141 work has been extended to the analysis of 
an induction motor with voltage control using thyristor- 
diode combination (thyrode) by Rahman and Shepherd [18] 
Arunachalam [191 presented the analysis of the same scheme 
m which state variable technique was used to obtain the 
initial vector and the final solution was obtained by the 
application of harmonic equivalent circuits Experimental 
observations of efficiency and power output of an induction 
motor controlled by delta connected thyristors at the open 
neutral of the stator winding was reported by Spooner [20] 
Hayashi [21] has performed the analysis of the same scheme 
using state variable approach 
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Paice [ 4 ] has established the fundamental laws relating 
the speed control of induction motors by voltage control 
schemes He also presented test results of eight different 
thyristor voltage control schemes to enable a comparative 
study A comparative study of circuits using six thyristors 
in different configurations has been reported by MdMurray l 22l 

142 Modelling of Induction Motors 

Two-axis models of induction motors described by Stanley 
equations 1 23] are commonly used for digital simulation [ 24,25] 
.Krause [261 has described the machine equations in an arbitrary 
reference frame One of the advantages of using the two-axis 
models is that the time dependence of the inductance coeffi- 
cients is eliminated Saito and Miyazawa [2/] described 
digital simulation of polyphase induction motors using a 
three phase model m which the machine equations are directly 
solved by finite difference methods for balanced and unbalanced 
conditions Using a three phase model, the digital simulation 
of an induction motor with non-sinusoidal excitation was 
presented by Sarkar and Berg 1281, The mam disadvantage of 
the three phase models [27,28] is that inverse of a time de- 
pendent matrix is to be computed at each step of integration 
Robertson and Hebbar [29] presented a digital model of an 
induction motor m which the machine behaviour is described 
directly in terms of the stator phase variables The variable 
stator-to-rotor mutual inductance is made constant by describing 
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the rotor variables m a stationary reference frame 

Ramshaw and Padiyar 1 303 described a generalized system 
model for all types of slip-ring machines including synchro- 
nous and induction motors A current source equivalent cir- 
cuit for an induction motor and its dynamical equations were 
presented by Padiyar and Ramshaw 1311 

14 3 Transient Analysis 

Numerous studies of the transient performance of indu- 
ction machines with constant applied voltage have been re- 
ported 132-361 and a wide range of performance aspects in- 
cluding starting, overspeeding, plugging and switching have 
been covered Performance of the induction motor with 
variable -frequency supply has been described m [ 37 1 Dyna- 
mic behaviour of the induction motor durin 0- balanced and un- 
balanced supply conditions were presented m 1381 using ana- 
log computers Transfer functions and dynamic response of 
induction motors with constant voltage have been presented 
139,401 on the basis of non-dimensional parameters 

Transient analysis of a voltage controlled induction 
motor drive with thyristor and thyrode controllers was pre- 
sented by Rahman and Shepherd ( 181 , using state space tech- 
niques The start up transient currents, torques and speed 
are obtained from digital simulation 

In the design of a speed controller for an induction 
motor with phase controlled tnacs, Kenly and Bose C 411 used 
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a very approximate first order model of an induction motor 
An experimental closed loop induction motor drive incorpora- 
ting phase controlled thyristors has been reported by 
Shepherd and Stanway [ 42 j 

1 5 OBJECTIVES AND SUMMARY OF THE VO BK 

The objectives of this thesis are to apply and develop 
new techniques for the analysis and simulation of the steady 
state and dynamic performance of voltage controlled induction 
motor drives Both open-loop and closed -loop operations for 
speed control are considered and analytical results are 
supplemented by experimental investigations A novel hybrid 
model of an induction motor is developed which is simple in 
structure, yet general enough to consider any set of stator 
terminal constraints introduced by thyristor switching This 
model is then used for the development of model for the dy- 
namic analysis of the drive 

A chapterwise summary of the work reported m this 
thesis is given below 

In Chapter 2, the steady state analysis of a voltage 
controlled 3-phase induction motor using back-to-back 
connected thyristors is described Boundary value approach 
143-46] is used for the analysis by which the steady state 
initial vector is obtained without an iterative procedure 
The machine equations are unaltered while forcing voltages 
are defined for each mode of operation Computed results are 


i 
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compared with experimental and other published results An 
equidistant pulse firing circuit [47] developed for the experi- 
mentation is briefly* discussed 

In Chapter 3, BVA is applied to the steady state analysis 
of delta controlled and thyrode controlled induction motor 
drives A comparative study of the performance of these 
drives is presented Digital simulation results 
showing the steady state performance are given 

In Chapter 4, two novel hybrid models are developed 
which are well suited for the steady state and transient 
analysis of induction motors under wide range of terminal 
constraints In the first model, rotor flux linkages on a 
synchronously rotating reference frame, are used and the 
stator is represented by a 3-phase circuit model with con- 
stant parameters The time varying mutual coupling between 
the stator and rotor windings is replaced by a dependent 
current source along the lines of 1 31 1 This enables the 
simulation of stator controlled drives by direct application 
of the network equations. In the second model, the stator 
variables are transformed into two phase variables on syn- 
chronously rotating reference frame and rotor variables are 
undisturbed so that any rotor terminal constraints can be 
directly incorporated Two examples are given to illustrate 
the applicability of the model for the digital simulation of 
thyristor controlled drives 



12 


Utilizing the circuit model developed in Chapter 4, a 
simplified analysis of steady state performance of a voltage 
controlled induction motor is described in Chapter 5* The 
analysis is based on the assumption that rotor flux linkages 
referred to a synchronously revolving frame are constants 
during steady state operation Analytical expressions for 
the stator current and voltage, during steady state operation, 
are derived Prom this, nonlinear algebraic equations for 
phase angle <f> and fundamental Pouner components of stator 
current are obtained as functions of triggering angle, supply 
voltage and rotor flux linkages By solving the set of iron- 
linear equations by Uev/ton-Raphson method, the steady state 
operating variables such as constant component of torque, 
rotor flux linkages, fundamental components of stator voltage 
and current and phase angle 4> are obtained for any given 
value of slip, triggering angle and supply voltage. The pre- 
dicted results match excellently -with the simulation results 
of Chapter 4. 

In Chapter 6, a third order linearized model of the drive 
about an opeiatmg point is developed to study the dynamic 
response due to perturbations in triggering angle, load 
torque and source voltage Stability of the drive based on 
eigenvalue analysis rs presented Utilizing the transfer 
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functions of the drive, a closed loop control scheme is 
designed and developed to regulate speed ki example illu- 
strating the design of the speed controller is given Ana- 
lytical and experimental results for the system response due 
to perburbat ions m load torque and reference voltage are 
furnished 

Chapter 7 summarizes the main contributions of the thesis 
and indicates the scope for further work m the area 



CHAPTER 2 


STEADY STATE ANALYSIS AND SIMULATION OF 3-PHASE INDUCTION 
MOTORS WITH VOLTAGE CONTROL BY SYMMETRICALLY TRIGGERED 

THYRISTORS 


2 1 INTRODUCTION 

One of the methods of speed control of induction motors is 
to vary the terminal voltage of the motor using thyristor swit- 
ching circuits £11J This type of speed control is inexpensive 
and reliable and is well suited for induction motors with fan- 
type loads or with high rotor resistance Different thyristor 
combinations can be suitably employed for the speed control of 
3-phase induction motors but a pair of back-to-back connected 
thyristors m each line is superior in many respects [43 » ! 

Although voltage control schemes are easy to realize m practice, 

i 

the analysis of such systems is considered to be difficult 

Steady-state analysis of a wye connected induction motor t 

l 

with voltage control using a pair of back-to-back connected , 

thyristors m each line is described in this ohapter. Boundary f 

i 

value approach (BVA) - ah, approach earlier found to be versatile I 
m relation + o the analysis of inverter-fed induction motors and [ 
h v d c systems [ 43 ^ 46 J „ is used for the analysis The 
steady state solution is obtained without altering the matrix > 
differential equation describing the motor for different modes 

\ 

of operation Change of mode only alters the excitation Open [ 
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circuit is simulated by applying a voltage, equal to the 
back-emf, across the open circuited phase The initial vector 
which will yield steady state solution is obtained without any 
iteration 

The basic assumptions made m the analysis are as follows 

a) The 3-phase voltage source is balanced and has zero source 
impedance 

b) All the thyristors are ideal 

c) All parameters of the machines are constant and magnetic 
saturation is neglected 

d) Rotor speed of the motor is constant during steady state 
operation 

e) MMP in the air-gap is assumed to be sinusoidal 

f) Hystens and eddy current losses are neglected 

2.2 DESCRIPTION OP THE SYSTEM 

Schematic diagram of an induction motor with voltage 
control is given in Pig 2 1 m which a pair of thyristors 
connected back-to-back is included in each line These six 
thyristors (T^-T^) are triggered symmetrically at intervals of 
it/3 to obtain a balanced voltage across the terminals of the 
machine By varying the triggering angle ( a ) (which is mea- 
sured from the zero crossing of the source voltage per phase), 
the effective voltage applied to the machine terminals can be 
varied to effect speed control 
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There are three distinct modes of operation for a vol- 
tage controlled wye connected induction motor -with isolated 
neutral These modes are identified m terms of the hold off 
angle y (see Pig 2 2), as follows 

1 ) Mode-1 0 < y < ir/3 

During this mode of operation either two or three thyri- 
stors will be conducting at any instant of time The sequence 
of thyristor conduction is T^ Tg, Tg Tg T.j , Tg T^ , Tg T^ Tg, 

mm rn m m mm m m m mm rp m m mm 
1 1 x 2 9 X 1 J 2 x 3 9 x 2 3’ x 2 3 X 4 J ‘ x 3 x 4* x 3 x 4 x 4 x 5* 

T 4 T 5 T 6> T 5 T 6 3,14 30 0,1 mrlae the 5 erioa when three thy - 

ristors are conducting all the three windings of the motor are 

connected to the power source while single phasing of the 

motor takes place when only two thyristors are conducting 

2) Mode-2 tt/3 < y < 2 w/3 

During this mode of operation, either one or two thyri- 

stors will be m conducting state When one thyristor alone 
is in conducting state, all phase currents will be identically 
zero because of the isolated neutral connection The sequence 
of thyristor conduction is Tg Tg, Tg, Tg T^ , T^ , T^ Tg, Tg, 

Tg Tg, Tg, Tg T^, T^ , T^ Tg, Tg, Tg Tg and so on 

3) Mode-3 2*/3 < Y < it 

During this mode of operation, all the stator currents 
will be identically zero and hence this mode is of no pra- 
ctical importance 
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2 3 MACHINE EQUATIONS 

The differential equations describing the 3-phase indu- 
ction motor expressed m per unit, -wherein the stator and 


rotor 

given 

variables are refeired to a stationary d-q 
m matrix form [2 3 , 48*49] by 

v ds 

| 

V X S D 0 
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V qs 
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v dr 
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Similar relationships hold good for the currents 
The matrix differential equation (21) is rewritten as 




v ds 
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where 


(2 7 ) 


V s1 = X s i ds + L dr 
V s2 * X s 1 qs + X m 1 qr 
V r1 = X m ^ds + X r x dr 
V r2 * 1 qs + X r 1 qr 
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The expression for the d-q currents m terms of the new variables 
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x 4. ’ < X r V S 1 ' X m V r1>/ (X 3 X r ' 

"qs - (X r V s2 ' \ - 4> (2f10) 

1 dr = < X s V r1 ' *m T al)/<*» x r ' X m> 

V (X s V r2 - *m 7 s2>/ (S 8 *r- 

Bauation ( 27 ) represents a linear, first order vector differen- 
tial equation of the form 

Bx = v -|A]x (2 11) 


where 


- ,v s1 V 7 r2> 


.1 


[V ds V qs 


0 0] 


Since the thyristors are triggered symmetrically , the voltage 
and current waveforms have half-wave and 3-phase symmetiy and 
hence there exists a boundary relationship between the varia- 
bles at u)t a 0 and wt = it/ 3 during steady state The boundary 
relationship as given m Appendix A is 


x 


j w t = 


it/ 3 ^ ~ j tut =0 


where 


T = 


1/2 

•*'3/2 

0 

0 


/ 3/2 

1/2 


0 


0 


0 

0 

1/2 

73/2 


(2 12 ) 


0 

0 

/ 3/2 

1/2 


* 



21 


The expression for the electromagnetic torque m per unit is 

T e * 2 / 3 V '"da V - 1 dr V ( 2 - 1 5) 

In the above expression for electromagnetic torque all 
the variables are to be m per unit The peak values 
of rated line to neutral voltage and rated line current 
are chosen as base quantities 

2 4 STEADY STATE SOLUTION BY BVA 

Equations (2,11) and (2 12) constitute a two-pomt 
boundary value problem for -which the steady state solution 
can be obtained as follows 

Referring to equation (2 11), one can rewrite this to include 
the unknown initial vector as 

Dx + Ax = v+ 5 (u>t ) x (214) 

where, 6 is the unit impulse function and the initial vector 
2* - 1 x 1i x 2j. x 3i x 4i> 

The equation (2 14) is split up into a number of component 
equations for convenience of numerical integration Super- 
position of the solution of these component equations, with 
all initial conditions set to zero, helps to evaluate the 
unknown initial vector x^ The component equations are * 
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Dx^+Ax^=v ( 215 a) 

D x 4 + A x 4 » 6(«t) e^ (2 15b) 

D x 2 + A x 2 = «(wt) e 2 (2 15c) 

D x^ + A x^ = <5 (art) (2 l5d) 

D x 4 + A x 4 = 6(wt) e 4 (2 15e) 

where 


e® * [1 0 0 0] , e g = ( 0 1 0 0 ] , e 5 = [ 0 0 1 0 ] , 

e^ = [0 0 0 1] and x°, x 1 , x 2 , ^ aad x 4 are 

components of x 

Lquations (2 15 a) to (2 15 e) are numerically solved using 
fourth order Runge-Kutta method over an interval of u/3, The 
origin is chosen as the current zero point of a-phase so that 
the foreing voltages for different operating modes can he 
easily defined Applying the principle of superposition, the 
final values of the solution at (wt «= tt/ 3) can he written as 

0 1 2 3 4 (1 

Xj? = Xf + ^“J 2 . ^ x 2l ^ Xj x Xjp ^ ^ ^ ' 

The subscript f stands for the final value at w t = Tr /3 

After rearranging the terms and making use of the boundary 
condition as given m equation (2 12), the initial vector 
becomes 

„-1 0 
— i * c 2f 


(2.17) 
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where 

[ CJ=s[T - X f J and t X f ] = [ xj ! I If t 1 

1 » f 

The initial vector obtained m (2 1 7 ) will lead to the steady 
state solution of (2 11) provided the forcing voltage vector 
v is defined for the period of integration 

2 5 DEFINITION" OF FORGING VOLTAGES 

Typical waveforms of d-q forcing voltages for mode-1 
and mode-2 operation of a voltage controlled motor during 
standstill, developed on the basis of equations (2 2 and 2 3) 
are given m Figs 2 3 and 2 4 For the solution of equation 
(2 15 a), we should know the voltages Vg g and v^ g for a period 
of ir/3 radians Starting from the current zero instant of 
a-phase the d-q voltages for various modes of operation are 
defined as follows 

2 5 1 Mode - 1 hold off angle (y) < ir/3 

During the interval 0 - tt/ 3, two distinct time periods 
exist , namely 

a) Time period I, 0 < cot v 

b) Time period II, y < wt < tr/3 

Time Period - I. 0 < mt < y 

During this period, current through a-phase is zero and a 
back-emf due to mutual coupling appears across a-phase winding 
Referring to Fig 2 5, the d-q voltages are 


& 

% 



X 73/2 


State of 
thyris- 
tor s 

f 


i 

1 0h 




ave forms during 







tages lor a 3-phase motor during 


v^ s = /*3/2 (back~emf) (2 18) 

v qs = / 3/2 V m cos(wt + *) (2 19) 

Time Period - II, Y < mt < tt/3 

All the phase windings are conducting during this period 
and the voltage across each phase winding will be same as the 
source voltage Hence, 

v dg = /~3/2"V m sin (wt + ♦) (2 20) 

v qs = ^ ^7 2 ooe(»t + 4>) (2 21) 

2 5 2 Mode - 2 y > ir/3 tut < 2ff/3 

Time period I, 0 < mt < (y - ff/3) 

Referring to Rig 2 6, the d-q voltages during this 
interval are given by 


v ds 

OK 

/3/2 (back-emf -j ) 


(2 22) 

v qs 


^3/2 (back-emfg) 


(2.23) 

Time period 

II, 

(y - it/ 3) < «t < 

*/3 


v ds 

3S 

J 3/2 (back-emf^ ) 


(2 24) 

v qs 

35 

^ V m cos ^ w ^ + <►) 


(2 25) 

It may 

be 

noted that the voltages 

defined 

above are not 


explicitly known since <$ and back-emf are still unknown 
These are obtained by an iterative procedure 
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2 6 STEADY STATE SIMULATION 

2 6 1 Algorithm for Simulation 

The step hy step procedure for the steady state analysis 
of a voltage controlled induction motor through digital 
simulation for a given slip and Y is described belo-w A 
simplified flow chart is shown in Eig 2 7 (a & b) 

St ep-1 A suitable value of <f> and back-emf are chosen 
A good approximation for <J> is the impedance angle of the 
passive equivalent circuit 

St ep-2 The mode of operation and accordingly the forcing 
voltages are identified With these values of forcing vol- 
tages and with zero initial vector, the system differential 
equation (2 15 a) is numerically integrated over one-sixth of a 
cycle In each step of integration during the period of 
current discontinuity, the bach-emf is suitably modified so 
that the resulting phase current is very near to zero It 
is assumed that the back-emf remains constant in each mte- 
\ gration step 

St ep-5 The solution of equations (2 15b to 2 I5e) are 
obtained numerically over one-sixth of a cycle 

St ep-4 The initial vector and the a-phase stator current 
for the initially assumed <j> are calculated using the equa- 
tions (2 17 & 2 10) 

( 

i St ep-5 If the initial value of a-phase current is non-zero, 

l 












P e t u r h 



Fig 2 7(b) How diagram nf the obrnutine 'hunge' - 
Sim jt jtion of a 4 phase induction motor 
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<f> is suitably modified and steps 2 and 4 are repeated 
until a-phase current becomes zero 

St ep-6 One more integration of equation (2 11) -with the 

computed initial vector produces the steady state solution 
of the variables over the specified period The phase varia- 
bles are calculfeed using the transformation eauations 

The value of 4> and initial value of back-emf are 
obtained by linear interpolation The initial vector and 
corresponding a-phase current are computed for two different 
values of $ and hence the correct value of <f> for which i-j a 
is zero is calculated by linear interpolation The initial 
guess value of 4 is taken as the impedance angle of the 
conventional passive equivalent circuit The initial value 
of back-emf is obtained m two steps In the first step, a 
back-eraf of 0 5 p u is assumed and resulting a-phase 
current is obtained after one integration step* In the 
second step, -0 5 p u voltage is chosen and current in a- 
phase is computed for the same interval Ey linear inter- 
polation, the back-emf corresponding to zero current is cal- 
culated and used for the computation 

During the period of current discont inuity (y)» the 
back-emf is to be modified, in each step of integration, 
such that the resulting phase current becomes oractically 
zero In the beginning of each step of integration, the 
back-emf is updated from its previous value by adding a small 
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increment and the integration is performed If the resultant 
phase current is not less than the preset accuracy level, 
the back-emf is further modified and integration is repeated 
The magnitude and sign of the phase current will give an indi- 
cation of the error m the chosen value of back-emf and 
accordingly the back-emf is moaified 

262 Time Varying Solution of Phase Variables Over One 
Complete Cycle 

Since the thyristors are triggered symmetrically , it is 
clear that the stator currents and voltages are similar but 
displaced by 120° The current and voltage waveforms also 
have half-wave symmetry Because of these properties, tho 
following relationships exist (see Appendix A) 

i 1a (wt + it/ 3) = - l^Cwt) (2 26) 

i 1a (u)t + 2 tt/ 3 ) = i 1o (»t) (2 27) 

i 1a (wt + it) = - i 1a Ut) (2 28) 

Similar relationships exist for the voltage waveforms. Using 

equations (2 26 to 2 2s), the current and voltage waveforms 
of a-phase can be constructed over one full cycle 

The instant aneous electromagnetic torque developed over 

I 

one-sixth of a cycle is obtained by substituting the time 
solution of d-q currents m equation (2 13) For the complete 

cycle, the waveform gets repeated six times The average 

value of the torque is obtained numerically 
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2 7 EXAMPLE 1 

In this example, a laboratory test machine whose para- 
meters are given in Table 2 1 is considered The steady 
state analysis is carried out by digital simulation outlined 
in the previous section and the results are experimentally 
v enf led 

271 Digital Simulation Results 

Computed waveforms of phase current, phase voltage and 
electromagnetic torque of the test machine over one cycle 
period are shown in Digs 2 8 and 2 9 Pig 2 8 pertains to 
mode-1 operation m which a slip = 0 106 and Y =35° are 
chosen for the computation Steady state waveforms during 
mode-2 operation are shown in Fig 2 9 for step = 1 and 
Y = 75° In Pig 2 10, the average value of electromagnetic 
torque is plotted against slip for constant values of y • 

272 Experimental Verification 
2721 Description 

Experiments are conducted on a 1 1 kW, 4 pole laboratory 
test machine subjected to voltage control A separately 
excited d c generator with a resistive load coupled to the 
induction motor is used for loading the drive Prom the 
output power of the d c generator, the torque developed by 
the induction motor is calculated, taking into account the 
losses m the d c machine Since the motor torque is a 
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Table 2 1 

DETAILS OF THE TEST MACHINES 

1 1 KW, 400 V, 50 Hz, 4 pole, 3 phase, wye-connect ed, squirrel 
case induction motor 

The measured parameters of the machine are 

R 1 « 9 15 Ohms J a 0 016 Rg~m 2 

= 3 268 Ohms D = 0 0053 N-m-sec/rad 

X sl = 8 526 Ohms 

X rl = 8 526 Ohms 

Xjjj = 1 26 432 Ohms 

R-j , R 2 , X s i» and Xjjj are the parameters of the passive 
equivalent circuit 

TABLE 2 2 

DETAILS OF THE MACHINE GIVEN IN REF [141 

1/3 hp, 220 V, 50 Hz, 4 pole wound-rotor machine 

The parameters of the machine expressed m per-unit using 
rated voltage as base volts and 375 W as base power are 

R 1 =0 0566 

R 2 = 0 1252 

X g = 1 0318 

X„ =1 0318 
r 

XL « o 9690 






current, voltage and 
j, V =3^* *=36.5^ 
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function of speed, the motor could he run stably at lov/er 
ranges of speeds where the torque-slip curve of the motor 
has negative slope The operation at these points were 
limited to short duration since the motor currents were exce- 
ssive The triggering angle a and hold off angle y are 

measured from the actual voltage and current waveforms A 

descript ion, of the triggering circuit for the control 
follows 

2 7 2 2 Triggering circuit 

An electronic triggering circuit (Fig 2 11) suitable 
for triggering six thyristors of a 3-phase voltage controlled 
system is developed The output waveforms at different 
stages are given in Fig 212 and a brief description of the 
circuit is given below 

The input signal to the circuit is obtained from one 
of the phases of the 3-phase source voltage An open loop 
comparator is used for obtaining the zero crossing of the 
input signal The square wave output of the zero-crossing 
detector is applied to a ramp generator The output of the 
ramp generator is compared with a reference voltage or an 
error signal for controlling the triggering angle a of the 
thyristors. The a can be adjusted over the full range by 
varying the reference voltage The pulse width of the com- 
porator output is different for different a’s With the 
help of a monostable , pulses of constant width for all 
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values of a are obtained By using a d c voltage blocking 
circuit, the d c component of monostable output is removed 
before it is applied to the phase looked loop (PLL) The PLL 
will provide necessary frequency multiplication and synchroni- 
zation The clock pulses (six times the input signal fre- 
quency) obtained from the PLL, are applied to a ring counter 
consisting of six I) flip-flops m cascade By means of a 
synchronizing pulse obtained from the monostable Mg (T'lg 
2 11 ), the ring counter is initially set to 100000 state 
There is provision for feedback input for automatic control 
of a The triggering angle correction, if any, is made m 
each cycle The circuit employs digital and linear [ 50,511 
integrated circuits The gate pulses are modulated and am- 
plified before it is applied to the thyristor gates through 
pulse transformers The gate current can be set between to 
10 mA and 500 mA to suit the gate characteristics of any 
type of thyristor The mam features of the triggering cir- 
cuit are 

1) suitable for 5-phase voltage controlled and bridge con- 
verter systems using thyristors 

2) self starting and self synchronizing capability 

5) suitable for various types of thyristors with different 
gate characteristics 

4) precise and wide range of & control 

iS 

5) inhibition of the gate pulse can be easily achieved 
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6) only single phase input is required 

7) due to the presence of synchronizing pulse throughout 
the operation, mal-operation of triggering sequence 
due to external electrical disturbances is corrected m 
every cycle 

8) minimum number of integrated circuit components are used 

2723 Experimental results 

Figure 2 10 shows, along with the computed results, the 
measured values of torque as a function of slip The actual 
voltage and current waveforms during mode-1 operation are 
given in Fig 2 13 When the triggering angle is ad ousted 
for mode-2 operation, the motor was practically stalling and 
the voltage and current waveforms during this mode of opera- 
tion are shown m Fig 2 14 

2 6 Example- 2 

This example is used to have a comparison of the results 
obtained by BVA with the results given m t 14 1 where state 
space technique is used for the analysis The average tor- 
que computed by BVA for various values of Y is given in 
Fig 2 15 The results from lipo 1 1 41 are also shown in the 
same figure The details of the machine are given in Table 
2 2 

2 9 Discussion 

The current, voltage and torque waveforms shown in Figs 
2 8 and 2 9 conform well with other published results in the 
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tip <? 1 3 Actual waveform o* voltage and current during 
nodc~1 



fig 2 14 Actual waveform of voltage and current during 
mod e«2 
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Pig 2 


Pig 2 

k 

\ 


s / \ 




13 Actual -waveform of voltage and current during 
mode-1 



14 Actual waveform of voltage and current during 
mode-2 


i 




literature The voltage -waveform is piecewise continuous 
and half-wave symmetric The electromagnetic torque contains 
a constant component with a predomment sixth harmonic The 
sixth harmonic torque is produced by the fifth and seventh 
harmonics m the stator current Harmonic contents of the 
stator current increase with increase m Y Nhen Y be- 
comes greater than it/ 3, the stator current becomes discon- 
tinuous m each half cycle with two current pulses Torque 
waveform also becomes discontinuous with six pulses in, each 
cycle The average toroue plotted against slip for constant 
values of Y (Pig 2 lo) has good agreement with the experi- 
mentally measured results The acrual voltage and current 
waveform (Pigs 2 13 and 2 14) also tally with the computed 
results 

Pigure 215 shows the torque speed characteristics com- 

\ 

puted by BVA for a motor given m (141 The results given in 
1 14) are also shown m Pig 2 15 to facilitate comparison 
It is observed that there is good agreement between the re- I 

suits obtained by BVA and those by state space technique 

Per unit quantities are used in the computation Peak \ 

t- 

value of the rated phase voltage and peak rated line current 
are taken as base quantities 

2 10 CONCLUSION 

Steady state analysis and simulation of a voltage controlled 
3-phase induction motor, have been presented using the boundary 
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value approach The main feature of this method is that the 
machine equations are unaltered for different terminal con- 
straints -while excitation functions are defined according to 
the mode of operation Hence, the mathematical formulation 
hecomes simple and general 


i 


# 1 * 



CHAPTER 3 


STEADY STATE ANALYSIS OE DELTA AND THYRODE CONTROLLED 
INDUCTION MOTORS BY BOUNDARY VALUE APPROACH 

3 1 INTRODUCTION 

Using BVA, the steady state analysis of a voltage contro- 
lled induction motor -with a pair of thyristors connected back- 
to-back in each phase -was described m the previous chapter. 

In order to demonstrate the generality of BVA, two more exam- 
ples of voltage control schemes are considered m thm chapter 
with, a) delta connected thyristors at the open star point of 
the stator winding, b) thyristor and diode (thyrode) combina- 
tion m the stator phases 

Although voltage control scheme with back-to-back conne- 
cted thyristor pairs is more efficient, other types of thyri- 
stor connections (4 } ^re sometimes used m practice when the 
prime consideration is not of efficiency but of economy* In 
a thyrode controlled scheme, three thyristors are replaced by 
three diodes (Eig 3 1) In a 'delta controller' three sy- 
mmetrically triggered, delta connected thyristors are inclu- 
ded at the open star point of the induction motor (Eig, 3*2) 
The mam disadvantage of delta and thyrode controllers over 
back-to-back connected thyristor pairs is the presence of 
predominant second harmonics m the stator currents These 
harmonics produce braking torque and also result in increased 
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iron and copper losses, which leads to a loss m the effi- 
ciency at low conduction angles 

3 2 S YMMET RIO ALLY TRIGGERED DELTA CONNECTED THYRISTOR 

CONTROLLER 

Three thyristors connected in delta are included at the 
open neutral of the stator winding The thyristors are tri- 
ggered symmetrically at intervals of 2 it/ 3 radians By con- 
trolling the triggering angle (a), the voltage applied to 
the motor terminals can he fully controlled Since the 
thyristors are connected m delta configuration, not more 
than two thyristors can conduct at any instant of time 
This is because of the fact that when two thyristors are 
conducting, the third one is reverse biased When two thyris- 
tors arc conducting, all the three stator windings carry current 
and voltage across each phase winding will be equal to the 
respective source voltage On the otherhand single phasing 
takes place when only one thyristor is conducting During 
this state of operation one of the stator rhases is open cir- 
cuited 

321 Modes of Operation 

Maximum period of conduction of each thyristor m a delta 
connected system is 4ir/3 radians If the delay angle or 
triggering angle (a) is zero, the thyristor R (Eig 3 2) 
starts conducting from a point 60° behind the positive zero 

1 

/ 

4 
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crossing of a-pnase voltage The conduction period (0) of 
each thyristor will he a function of a, slip and the machine 
parameters If the value of $ ranges between 2 ir/3 and 4 ir/3, 
either three or two windings of the induction motor are conne- 
cted to the power source This mode of operation is denoted as 
either 2/1 mode or mode-1 operation The conduction period of 
each thyristor and typical current and voltage waveforms are 
shown m Pig 3 3 for 2/1 mode of operation m which <* , Y 
and $ are indicated The thyristors R,Y and B are triggered 
at a, 2*/3 + a and 4^/3 + a respectively The conducting 
thyristors go into the blocking state at the natural current 
zero The sequence of thyristor conduction m this mode is 
B-R, R, R-Y, Y, Y-B, B, B-R and so on 


\Jhen the value of B lies between 0 and 2 it/3 either one 
thyristor or none will be m conducting state Single phasing 
takes place when only one thyristor is conducting When none 
of the thyristors is conducting, all the stator currents will 
be identically zero but voltages may be present across the 
windings (induced by the rotor currents) This mode of opera- 
tion is denoted as either 1/0 mode or mode-2 operation Ty- 
pical voltage and current waveforms are given m Pig 3.4 
Two current pulses per cycle appear m each phase winding 
and this mode of operation is not having any practical im- 
portance because the electromagnetic torque produced will be 

***• A 85 ® 

'-”*****&* * 


quite small 
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5 2 2 Symmetry of the System 

Since the three thyristors are symmetrically triggered 
at an interval of 2 ir/ 3 , the currents and voltages have 3 -phase 
symmetry However, unlike the case of six thyristor conne- 
ction, here half wave symmetry does not exist Because of 
the 3 -phase symmetry, hold off angle Y , phase angle <t> , delay 
angle a associated with the turn off of each thyristor are 
identical Accordingly a solution for any 2tr/3 interval is 
sufficient to uniquely define all steady state variables 
over the entire period A boundary relationship between the 
steady state variables at an interval of 2 ir /3 is given m 
Appendix B 


3 2 3 Machine .Equations 


The equations (2 1 ) to (2 17) derived m Ohapter 2 are 
directly applicable to the present case except (2 12) de- 
fining the boundary condition In the present case a boun- 
dary relationship exists during steady state for an interval 
of 2 f 3 radians The boundary relationship (see Appendix B) 
is given by 


where 


T. 



= 2 tt /3 - 

T 1 -|wt 

* 0 



- 1/2 

✓ 3/2 

0 

0 


- / 3/2 

- 1/2 

0 

0 


0 

0 

- 1/2 

/ 3/2 


0 

0 

-/ 3/2 

- 1/2 


(3 1 ) 
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The component equations of the machine differential 

equation (2 15 a - 2 I 5 e) are solved numerically over a 

period of 2 ir /3 The forcing function needed for the solution 

of (2 15a) are described m Section 3 2 4 The final value 

of the variables over an interval of 2 rr /3 is 

0 1 2 3 ,4 / -z o\ 

if = Sf + x 1i if + x 2i if + x 5f if + x 4f if (3 2) 

Combining equations(3 1) and ( 32 ), the initial vector can 
be obtained as 



■where 

Cg * T^ - and is a 4x4 matny given by 

_ r 1 2 ' 3 1 4, 

4f B 1 , if ' Xf 1 Xf 1 


One more integration of the system differential equation 
(2 15 a) over an interval of 2*nr/ 3 with the initial vector as 
computed m (3 3 ) provides the steady state solution 

324 Definition of forcing Voltages 

forcing voltages v ^ 0 and v qg are to be defined over a 
period of 2 n /3 for the solution of system differential equa- 
tion (2 15 a) for the purpose of analysis, the reference 
point for defining the forcing voltages is chosen as the 
instant at which the thyristor R extinguishes (figs 3 3 and 
3 4) As explained m Section 321 , there are two modes of 
operation depending upon the conduction period of the thyristors 
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3 2 4 1 2/1 mode of operation 

The d-q voltages as defined in equations (2 2 and 2 3) 
are plotted in Fig 3 3 for 2/1 mode of operation There are 
t-wo distinct time periods during the interval 2it/ 3 

a) Time period I, 0 < cot < y 

b) Time period II, Y < cot < 2 it/ 3 

Time Period I. 0 < cot < y 

At the beginning of the interval the thyristor R goes to 
blocking mode resulting in zero current in phase-a A back- 
emf of unknown magnitude appears across the winding* Thyri- 
stor I continues to conduct resulting in single phasing of 
the motor Referring to Fig 3 3, the d-q voltages are 

v dB * ^3/ 2 (back-emf) (3 4) 

v qs = % cos (u»t +4) (3 5) 

Time Period II, Y < mt < 2ir/3 

At ut = y the thyristor B is turned on resulting in 
3 phase conduction The voltage across each stator winding 
will be same as the respective source voltages The d-q vol- 
tages are 

v djg = - \ sin (cot + ♦) 

v qg = - ^3/2 cos (cot + <J>) 


£ 


(3 6 ) 
(3 7) 
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3 2 4 2 1/0 mode of operation 

The conduction period of each thyristor during this mode 
of operation is less than 2 v/3 radians The machine -will he 
subjected to single phasing and off state alternately Two 
time periods can he identified (fig 3 4) for the purpose of 
defining the forcing voltages 

a) Time Period I, 0 < wt < (y - 2ir/3) 

At the beginning of this interval, the thyristor R goes 
to off state Currents m all the phases will be zero as no 
thyristor is conducting Back-emf may appear across the 
phase windings due to the coupling with the rotor currents 
The forcing voltages during this period is 


v ds = /3 / 2 

(back-emf .j ) 

(3 8) 

v qs 

(back-emf 2 ) 

(3 9) 


b) Time Period II, (y - 2ir/3) < mt < 2*/3 

At ut = y , thyristor I is turned on resulting in 
currents m b and c phases The d~q voltages are (fig 3 4) 

v ds “ (back-emf^) (3 10) 

V = -*^3/2 v, cos (wt + <j>) (3 11) 

qs in 

3 2 5 Steady State Simulation 

After having defined the forcing voltages over a period 
of 2tr/3, the solution of system differential equations 


d 

j 
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(2 15 a - 2 , 15 e) can be obtained numerically for a given slip 
and Y as described xn Section 2 6 of Chapter 2 The final 
values of the solution obtained for a period of 2ir/3 are sub- 
stituted m equation (3 3) to obtain the initial vector This 
initial vector will directly lead to the steady state solu- 
tion 

The value of <t> is obtained by linear interpolat ion 
technique The back-emf in each step of integration during 
the period of Y is computed as described in Section 2 6 of 
Chapter 2 

326 Solution of Phase Variables Over One Cycle 

Prom the computed results of the state variables V -j , 

V g 2 > » ^r2* P^ ase currents i*j a > x^ and i^ c can be cal- 
culated using equations (2 10) and (24-26) Because of the 
3-phnse symmetry of the phase variables, the following rela- 
tionships can be easily derived (see Appendix B) 

1 1a^ <1> ^ + = Hc (ut) (3.12) 

i 1a (<*t + 4V3) = i J j- b (wt ) (3 13) 

Similar relationships exist for the stator voltage waveforms 
The electromagnetic torque gets repeated three times over 
one cycle 

327 Example 

Steady state waveforms of stator current stator 

voltage v and electromagnetic torque T are computed with 
as © 
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the forcing voltages defined in. Section 3 2,4 The algorithm 
given m Section 2 6 of Chapter 2 is used for the simulation 
Pig 3 5 shows the computed steady state waveforms of current, 
voltage and torque for a hold off angle ( y ) of 50° and slip 
of 0 15 The computed results for y = 100° and slip =02 
are shown m Pig 3 6 The dissimilar pattern of voltage 
and current waveforms during negative and positive half -cycle 
is evident from these results The torque waveform has a 
constant term with a predominant third harmonic The magni- 
tude of the harmonic component changes with y and this is 
evident from the torque waveforms given in Pigs 3.5 and 3 6 
It can he seen from Pig 3 5 that the peak: value of current 
reaches about 4 p u at a slip of 0 15 This can cause ex- 
cessive heating and damage the windings if the motor is run 
for longer duration at this value of slip The average 
value of the developed torque for constant values of y is 
shown m Pig 3 7 The parameters of the motor are given m 
Table 2 1 of Chapter 2 

3 3 SYMMETRICALLY TRIGGERED THYRODE CONTROLLER 

Schematic diagram of a voltage controlled 3-phase indu- 
ction motor with back-to-back connected thyfistar-diode combi- 
nation m each phase is shown in Pig 3 1 The three thyristors 
T^ , Tg and T^ are triggered symmetrically at intervals of 
2^/3 so as to obtain a balanced 3-phase voltage across the 
machine terminals The input voltage to the machine terminals 
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can "be controlled 1 ay varying the triggering angle a Because 
of the symmetrical triggering of the thyristors, the machine 
voltages and currents have 3-phase symmetry hut there -will not 
he any half wave symmetry due to the presence of diodes m 
each branch Depending upon the hold off angle Y of phase 
current, the motor operation is classified into two distinct 
modes If Y is less than 2 tt/ 3, the machine is subjected to 
three phase and single phase operation alternately Three 
devices - either two thyristors and one diode or two diodes 
and one thyristor - will le conducting when 3-phase opera- 
tion takes place During single phasing only t^o devices - 
a thyristor and diode - conduct Accordingly this mode of 
operation is identified as 3/2 mode When the value of y is 
greater than 2 tt/ 3 , the machine is subjected to single phasing 
and open circuiting alternately This mode of operation is 
denoted as 2/0 mode The method of steady state solution 
of motor variables for both modes of operation are exactly 
identical to that of the delta controller described in 
Section 3 1 Equations (3 1 - 3 3) are also applicable to 

the present case 

3 3 1 Definition of Forcing Voltages 

For the purpose of defining the forcing voltages, the 
current zero instant of a-phase is chosen as reference point 
The d-g voltages for both modes of operation are to be defined 
over a period of 2ir/3 to obtain the solution of system differe- 
tial equation (2 15 a) 
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Mode 2/3 0 < y < 2w/3 

Referring to Pig 3 8, the d-q voltages can be 
defined as 

a) Time period I, 0 < <ot < y 

V ,= back-emf 
ds 

v qg - * / 3/2 cos (tot + 4>) 

b) Time period II, y < wt < 2 tt/3 

v ds “ V m sin ^ t + ♦) 

v qg * ^3/2 V m cos (tot + (f>) 

Mode 2/0 Y > 2ir/3 

The d-q voltages as given m Pig 3 9 are 

a) Time period I, 0 K tot < (y - 2ir/3) 

v^ s » ^3/2 (back-emf^ ) 

v qg =s ^3/2 (back-emf 2 ) 

b) Time period II, (Y - 2 tt/ 3) < tot < 2ir/3 

v^ g a ^3/2 (back-emf) 
v qg - *^3/2 V m cos( w t + 4) 

3 3 2 Steady State Simulation 



readily 

(3 14) 
(3 15) 

(3*16) 
(3 17) 


(3 18) 
(3,1 9) 

(3 20 ) 
( 3 . 21 ) 


The procedure to obtain the steady state waveforms of 
current, voltage and torque are exactly similar to that of a 
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delta controller described m Section 3 1 The forcing 
voltages needed for the solution of the system differential 
equation (2 15a) are given m equations (3 14 - 3 21) 

3 3 3 Example 

Steady state waveforms of phase current i^ a phase vol- 
tage v and electromagnetic torque obtained from the com- 
puter simulation are shown m Pig 3 10, for 2/3 mode of 
operation (Y = 40°, slip = 0 1) It can be observed from 

Pig 3 10 that there is no half-wave symmetry for the current 
and voltage waveforms The to roue waveform has a predominant 
third harmonic component as m the case of the delta contro- 
ller. Pig 3 11 shows the computed torque slip characteri- 
stics for constant values of y 

3 4 DISCUSSION 

Prom the computed results, it can be seen that the tor- 
que in a delta and thyrode controlled machines has a pre- 
dominant third harmonic component The third harmonic torque 
is produced by the second and fourth harmonic components of 
the stator current Although the shape of the current and 
the voltage waveforms (Pigs 3 5 and 3 10 ) of the delta and 
thyrode controllers appears to be different, the r*m,s. 
current and average torque for a given slip and Y are exactly 
same for both cases This fact can be observed from Table 3 1 
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Fig 3 10 Computed results of a thyrode 
controlled motor (slip= 1 > i « 4C?> 
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TABLE 3 1 

STEADY STATE SIMULATION RESULTS OE THYRODE AND 

DELTA CONTROLLERS 


SI 

No 

y 

slip 


Thy rode 

Delta 


m s 

T 

ave 

I 

r m s 

T 

ave 

1 

10 

0 05 

1 587 

1 483 

1 587 

1 483 

2 

10 

0 10 

2 387 

2 032 

2 387 

2 032 

3 

10 

0 20 

3 371 

2 140 

3 371 

2 140 

4 

10 

0 30 

3 858 

1 900 

3 858 

1 900 

5 

10 

0 40 

4 140 

1 640 

4 140 

1 640 

6 

30 

0 05 

1 550 

1 310 

1 550 

1 310 

7 

30 

0 10 

2 250 

1 780 

2 250 

1 780 

8 

30 

0 20 

3 090 

1 760 

3 090 

1 760 

9 

30 

0 30 

3 490 

1 490 

3.490 

1 490 

10 

30 

0 40 

3 715 

1 260 

3 715 

1 260 

11 

50 

0 05 

1 467 

1 100 

1 467 

1 10 

12 

50 

0 10 

2 110 

1 449 

2 1 1 0 

1 449 

13 

50 

0 20 

2 838 

1 338 

2.838 

1 338 

14 

50 

0 30 

3 091 

1 086 

3.091 

1 ,086 

15 

50 

0 40 

3 227 

0 893 

3 227 

0 893 

16 

70 

0 05 

1 376 

0 917 

1.376 

0 917 

17 

70 

Q f 10 

1 997 

1 087 

1 997 

1,087 

18 

70 

0 20 

2 510 

0 937 

2 5lO 

0 937 

19 

70 

0 30 

2 694 

0 743 

2.694 

0 743 

20 

70 

0 40 

2 790 

0 600 

2 790 

0 600 


C cor\td, ) 
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SI 

No 

y 

slip 

Thy rode 

Delta 

m s 

T 

ave 

I 

r m s 

T 

ave 

21 

90 

0 05 

1 379 

0 724 

1 379 

0 724 

22 

90 

0 10 

1 817 

0 790 

1 817 

0 790 

23 

90 

0 20 

2 280 

0 680 

2 280 

0 680 

24 

90 

0 30 

2 359 

0 507 

2 359 

0 507 

25 

90 

0.40 

2,386 

0*392 

2 386 

0.392 

26 

110 

0 05 

1 354 

0 550 

1 354 

0.550 

27 

110 

0 10 

1 696 

0 649 

1 696 

0 649 

28 

110 

0 20 

1 978 

0 508 

1 978 

0 508 

29 

110 

0 30 

2 060 

0 351 

2 060 

1 351 

30 

110 

0 40 

2 130 

0 284 

2 130 

0 284 


The current and torque are m per unit 
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where the computed results are given for both controllers 
The torque-speed characteristics (Pigs 3 7 and 3 11) for 
both controllers are also identical The forcing voltages 
defined m Sections 3 2 4 and 3 3 1 also indicate that the 
forcing voltages m both cases are the same except for a 
phase shift of 180° The motor is subjected to same dura- 
tion of single-phase and three-phase operations for a given 
Y and slip and hence, both controllers should result m the 
same motor performance 

Although there is some saving m terns of the thyristor 
cost, the current rating of the thyristors is higher than 
that of six thyristor scheme In addition, the motor is 
also subjected to higher heating in delta and thyrode con- 
trollers The ratio of the motor current (r m s ) of the 
delta (or thyrode) controller to that of a six thyristor con- 
troller shown m Pig 3 12 clearly indicates that the current 
ratio increases with hold off angle Y The reason for the 
larger values of current in delta and thyrode controller is 
the presence of even harmonics m the motor current 

3 5 CONCLUSION 

BVA Is extended to the analysis of two more cases of 
voltage controlled 3-phase induction motors m which 
different thyristor combinations are used In the first 
case, symmetrically triggered delta connected thyristors 
are used for the voltage control In the second case 
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thyristor-diode combination is used for the control In 
both cases the effective value of the motor voltage can be 
controlled over full range by varying the triggering angle 
Because of the symmetry of the system, a solution for the 
time period of a third of a cycle completely defines the 
motor operation It is observed from the computed results 
that performances of 3-phase induction motor controlled by 
symmetrically triggered delta connected thyristors and thyn- 
stor-diode combination are identical. 
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CHAPTER 4 


DIGITAL SIMULATION" OE THREE PRASE INDUCTION" MOTORS 
USING HYBRID MODELS 


4 1 INTRODUCTION 

In this chapter two models of an induction motor, that 
are well suited for the transient and steady state analysis 
of three phase induction motors, under a wide range of ter- 
minal constraints, are given The ma^or feature, of the models 
presented, is the development of an equivalent circuit for 
the stator or rotor windings, m which a dependent current 
source represeats the coupling "between the stator and the 
rotor This facilitates the formulation of the state equa- 
tions for the motor with stator or rotor connected to any 
arbitrary network 

One of the models presented is particularly suitable 
for stator controlled motors while the other is well suited 
for rotor controlled ones In the first model, the 3-phase 
rotor variables are transformed into equivalent 2-phase 
variables while the 3-phase stator variables are undisturbed 
This model gives the best approach f 5 2| to the class of 
problems m which mixed terminal conditions exist In the 
second model, the stator variables are transformed into 
equivalent 2-phase variables and rotor variables are undis- 
turbed so that any rotor terminal constraints can be directly 


? % 
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incorporated in the model Two examples are provided to 
illustrate the applicability of one of the models for wide range 
of terminal constraints, 

4 2 DERIVATION OP THE MODELS 

The following assumptions are made in setting up the 
machine equations 

a) no magnetic saturation 

b) uniform air-gap 

c) negligible hysteresis and eddy current losses 

d) sinusoidal distribution of air-gap flux 

4 2 1 Pirst Model With Equivalent Circuit for the Stator 

The rotor is assumed to have two short circuited windings 
on a - 8 axes rotating at a speed of u r Referring to Pig 
4 1 the stator and rotor voltages can be written as 


Z.* * io + P 1 


V as R 1 + T) l!) 

— r r — r v 


(4 1) 

( 4 . 2 ) 


where 


r 


R. 


R 1 0 0 

0 R 1 0 

0 0 R, 


, R r = 


r 2 0 

0 R, 


T 

* V o 

9 — S 


«* Vj 


as bs cs 


1, 


v r - I 0 01, i g = l i-| a 1 i c I » * * x 2a 1 2$ ^ 


Is -<*1a *1b tx 


i.T 
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The stator and rotor flux linkages are gi/en hy 

* L i e + M 1 

s —s — r 

= I i + i 

~r r — r + is 

where 



1/^ 

P 1 Pi 


Ji z 0 

P - 

Pi 

Li L 1t 

, L r - 

0 Lo 

i 

Pi 

Pi P . 


L <i 

and 


cos 9 

1 r 

sm 

e r 

M * ] 

p2 

i 

cos (0 r -.2Tr/3) sm(9 r -2ir/3) 



cos (0 r +2ir/3) sm(0 r +2n/3) 

-Hr* 


Prom equation (4 4), l can be expressed as 


-1 _i T 

l = {1 1 ty - L 1 M l > 
— r r — r r is 


Substituting equation (4*5) in (4 3)* we get 
i|> * (1 - M I" 1 M T ) i + M L" 1 ♦. 

— S R T ' — R T* — . 


Let 


1" . (I B - M IT 1 M T ) 


i - l"' 1 m x,; 1 *_ 

—s s r r r 


r -r 


(4 3) 
(4 4) 


<4 5) 


(4.6) 


(4 7) 
(4 8) 


* 
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Then eqn (4 6) can he written as 

= L o Ci c + i e ) 

—s s —s — s 

Equation (4 7) can he expanded as 



1" 

aa 

=5b 


K - 

S 

*ib 


L ab 


L ab 

“Sb 

L aa 


where 


L" 

aa 


L 1 " L 12/ L 2’ L ab * L 11 + l 12/^ 2L 2^ 


On expanding (4 

8), 




cos 0 r 

1 

sm e r 

Is - 1 - 

cos(6 r -2tr/3) 

■ + 1 $ 

sm(e r -2ir/3) 


cos(9 r +2iy , 3) 


sin (6 r +2V3) 


where 


L 12 ^ , *12 


(4 9) 


(4 10 ) 


(4 11) 


13 = 


I" - L" 

3,9, 3D 


I a M L 2 Tf I $ c Lpq r ^2B' L 1 

The stator voltage equation (4 1) can he rewritten as 

Z. - * B ±8 *<i. - J.) (4<12) 


Equation (4 12) can he represented by means of an equivalent 

circuit (Pig 4 2), which will he useful for the formulation 
of the state equations for any stator terminal constraints 


* 
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Comb ii mg eqns (4 2) and (4 5), we get 
p _$ ss -R (IT 1 - IT 1 M T i ) 

* -ip 3? T X* 


(4 13) 


On expanding, 


*2 a 

r 2 

*2ci 

R 2^12 ; 

cos0 r cos(e r -2Tr/3) cos(0 r +2*/3) 

_+29 j 


,*2bJ 

+ V 

sm0 sm(e -2 tt/ 3) sm( 0 + 2ir/3) 

J. jL -J- — 


1 1 a 
1 1b 


(4 14) 


While equation (4 14) can be directly used for digital simu- 
lation, it will be advantageous if the rotor flux linkages m 
eons (4 11) and (4 14) are transformed into d-q variables 
(referred to a synchronously rotating reference frame) so that 
the transformed variables become constants during steady state 
operation with sinusoidal excitation Even for non-sinu'-’oidal 
excitation such as due to the operation with periodic switching 
of supply, it can be shown that the rotor flux linkages 
and if> 2 q ( in the synchronously revolving frame) are practically 
constants in steady state (at constant speed) operation 


The transf ormat ion to d-q axis is defined by 


♦ 

2a 


cos(e o - e r ) 

sin(0 o -0 r ) 


* 24 

^2 0 


-sm(e 0 - e r ) 

cos (e 0 -e r ) 


" , 2q 


(4 15) 
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In teems of the d-q variables, equation (4 11) "becomes 



cos 0 


U*. 

sin ® 


0 


0 

= h 

Cos( 9 q — 2 tt/ 3 ) 

+ 1 2 

sin(0 o -2ir/3) 


cos(© 0 +2it/3) 


sin(0 o +2iT/3) 


(4 16) 


■where 


llA r> 

-h “ ijnf ^2d ’ 


T 12 * 

I 2 * v 2q 


Using (4 15), eqn (4 14) can "be transformed as 
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-R 2 /L 2 -Sw 
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1 1d 
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S u> 
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(4.17) 


where 


i — -» 

1 1d 


•1q 


* / 2/3 


cos 8 q cos(6 0 -2it/3) cos(0 o +2tr/3) *j 

sin 0 Q sin(6 o -2ir/3) sin(0 o +2ir/3) 


X 1 a 
x 1b 

"He 


(4*18) 


Equation (4 17) represents the rotor flux linkage equation 
for any type of stator terminal conditions The instantaneous 
torque (as derived m Appendix 0) and speed equations are given 

ty 

L. 

J 2 

and 

P s . . ( *« : T V. + H4L£ (4 so) 


^(1 14 *2, - hq 

(3?- “ Tl 


(4 19) 
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Equations (4 12 ), (4 17), (4 19) and (4 20) completely des- 
cribe the induction machine under dynamic conditions 

The equivalent circuit (Pig 4 . 2 ) representing a general 
induction machine can be further simplified if the stator 
■windings of the motor are wye connected with isolated neutral 
Since the sum of the stator currents and their derivatives 
are zero at any instant of time and by virtue of eqn (4 16), 
the stator voltage eqn (4 12) can be rewritten as 


r n 

v «is 


r -i 
Ha 


** 

x 1 a 

4 

I 1a 

v bs 

- R 1 ’ 

1 1 b 

+ L^’ p 

Hb 

4 

hb 

V d S 


_H<v 


x 1 c 

4 

ha 


The equation (4 21 ) is represented by means of an equivalent 
circuit (Pig 4 3) m which there is no mutual coupling bet- 
ween the stator windings In Appendix D, the relationship 
between the parameters of conventional passive equivalent 
circuit and those of the active equivalent circuits (Pigs 
4 2 and 4 3) are given 

4 2 2 Second Model with Equivalent Circuit for the Rotor 

In this model, the rotor variables are undisturbed while 
the stator variables, are first transformed into a- $ varia- 
bles Referring to Pig 4.4, the machine equations are 


V* 

— s 

» R r 
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— r 

= R r 
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— r 
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(4 22 ) 

(4 23) 
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Rg 4 3 Circuit model of a wye- 
connected stator 


Fig 4 4 Schematic diagram of 
an induction motor 
with 2 phase stator 
and 3 phase rotor 
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■where 
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The relationship between 3-phase voltages and a 
given by 
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(4 24) 
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From equation (4 25), the stator current i^ can he expressed 


li = { Li* 1 jfe* - L' “ 1 II' T 1« } 


(4 27) 


Eliminating i' from equation (4 26) using (4 27) 

-~s 
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(4 30) 


The rotor voltage eauation (4 23) can he rewritten as 
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(4 31) 


The stator flux linkage equation m terms of the rotor currents 


i» can he obtained from equations (4 22) and (4 27 ) as 
— r 


p ^ ’ a v 1 - R’ L'" 1 V' + R1 l'” 1 M' T ii 
p -s — s s s — s s s -r 


(4 32) 


Aftei transforming the ct 
using the relationship 


6 components to d-Q components 
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(4 34) 
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Por a wye --connect ed rotor winding with isolated neutral, 
the rotor voltage equation (4 31) becomes 


~ V 2a~ 


i 2a 


1 2 a * I 2a 

& 

C\J 

P> 

- R Z 

x 2b 

+ I>2 P 

1 2b + I 2b 

V 2c 


_ l2c _ 


1 2c + I 2 C 

— 


(4 36) 


equation (4 36) is represented by an equivalent circuit 
(Pig 4 5) where Ig , and l 2 c are ^ e P en dent current sources 

which are functions of stator flux linkages ana The 

instant nneous torque (as derived m Appendix E) and the speed 
equations are 


T = /vT 


p8 = - 


(T. 


tip (l 2q " X 2d q ) 

i ~ V (1-S) D 
wj + J 


(4 37) 
(4 38) 


4 3 EXAMPLES OP STATOR CONTROLLED MOTORS 

4 3 1 Example 1 Digital Simulation of a Wye-Connected 
Induction Motor 

The schematic diagram* shown in Pig 2 1 (Chapter 2) is 
considered for digital simulation The terminal voltage of 
the motor is varied by symmetrically triggered, back-to-back 
connected thyristor pairs m each line Por dynamic and 
steady state analysis of such a system, the different possible 
modes of operation are 
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1) the 3-phase source is directly connected to the motor 
terminals allowing currents m all the three lines 

I 

2) line a' a is open while other two lines carry currents 

3) line b'b is open while lines a’ a and c’o carry currents 

4) line c’c is open while lines a' a and h'h carry currents 

5) currents m all the three line® are zero 

This mode of operation will occur only when not more 

! 

than one thyristor is in conducting mode I 

t 

4 3 11 Formulation j 

The objective is to obtain a set of equations corres- 

ponding to each mode of operation from network analysis j 

\ 

For a particular case when the supply is assumed to be a 
voltage source m series with an inductance, the network on , 

the stator side is shown m Fig 4 6(a) Here it is assumed 
that the thyristors m all phases are conducting 

{ 

There are three possible meshes for the network shown 
m Fig 4 6(a) although for a given mode, not more than two 
mesh currents are independent The equations for each mode 
would be given m terms of mesh currents For each of the 
five modes of operation, the state equation can be formulated 
as follows 

Mode-1 

let and i 2 be the independent mesh currents The 
mesh equations m matrix form are given by 

j 
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a'n~ b'n 2(l g +L^) ~(l g +lp 


V -V 

b'n c'n 


(VI?) 2 ( 1 S+ I|) P [i 2 J + 

fa -*;1 f ',.1 


L!| 2 1»Y 


2R^ — R^ i “1 


-R 1 2R 1 


(4 39) 


Equation (4 39) can be simplified to 


i-i i-i 10 I 4 

<V L 1> ? -J4 -1? P 1a 

■ O 1 O * I -1- J 


(4 40) 


Equations (4.39) and (4 40) are obtained using the relations 

v a’n + v b'n + v c'n * 0 aild X 1 a + I 1'b + T 1c * 0 The relat;L 
ship between the phase currents and mesh currents are 


1 0 0 

-1 1 0 

0-10 


(4 41) 


where is the third mesh current which is not considered 
m this mode of operation 

Mode-2 

When the stator current i 1a becomes zero, the thyristor 
in a-phase enters the blocking state During this mode, back- 
emf will appear across the a-phase winding due to induced 
rotor currents From the equivalent circuit (Eig 4.6(h)), 
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the expression for the back-emf is 


v a(back-emf ) ^ 


(4 42) 


and the state equation is 


+LL' )' + v b’n " v c’n ^ ^ 

S 1 


The phase currents are 


0 0 0 


0 1 0 


o-l 0 


Mode-3 


(4 44) 


Prom the equivalent circuit (Pig 4 6(c)), the state 
equation for this mode of operation is 


J+JEJ7 ^“ 2R 1 1 3"' I '1 p ^ ;i '1a"' I 1c^ + v a'n * v c'n ^ ^ 


The stator phase currents are 


0 0 0 


(4 46) 


0 0-1 


Mode-4 


This mode of operation is identical to those of 2 and 3 
except that the phase c is open-circuited, The state equation 
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and the stator currents are (ref Fig 4 6d) 
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Mode-5 


(4 48) 


(4 49) 


None of the stator windings carry current during this 
mode of operation hut back-emf may appear across the stator 
windings duo to the presence of rotor currents Referring 


to Fig 

4 6e 






i H * l. 

1 a 

1b * X 1 

c 

0 


(4 

50) 

v a (back- 

- emf ) 

hl{ 

P 

U a 



v b (back- 

-emf ) 


P 

I 1T) 

(4 

51) 

v c(back 

-emf) 

1 " 

P 

I 1c 




4 3 12 Simulation and results 

The rotor flux linkages ^ 2d and ^ 2q , slip S and the mesh 
currents , i 2 and i^ are chosen as state variables The 
state equation defining <f> 2d , ^ 2 q and S are derived m equa- 
tions (4 17) and (4 20) The state equations of mesh currents 
consist of terms with derivatives of dependent current 
sources which can be obtained from equation (4 16) as 
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Ii 

P x 1 a = Xptr {cos 0 O p t^+sm 9 o P *2q + °^ C0S 6 0 * 2q“ sin ®o W } 

(4.52) 

Xi 

pl lb « {cos(0 o -2rr/3) p t|> 2d + sm(8 0 -2ir/3) p t|> 2 q + 

u(cos(0 o -2ir/3) ij> 2g - sin(e o -2ir/5) <!» 2d )} U 53) 

and 

L 

pi. = y~£r ^os(e +2V/3) p t|»2d + sin (° 0 +27r / 5 ) P ^2q + 

w(cos(9 o +2it/3) - sm( 0 o +2V3) ^ 2q ) * 54) 

The state equations for the stator -which depend upon the 
particular mode of operation, are solved simultaneously along 
with the rotor flux linkage equations (4 17) and speed equa- 
tion (4 20) on a digital computer using fourth order Runge- 
Kutta method for a chosen triggering angle t* and load torque 

T A sten size of 0 00025 sec is chosen for computation 
m v 

Linear interpolation is used to obtain instants for the cu- 
rrent zero and its detection of change of mode initiated by 
the blocking of a thyristor Results of the simulation are 
shown in Figs 4 7 to 4 11 Transient starting torque for J 

a = 80° with zero load torque is shown in Fig 4 7 Figure 
4 8 shows the starting characteristics (speed versus time) 
for different triggering angles where the high frequency 
oscillations at the starting period are not shown It can be 
observed from Figs 4 7 and 4 8 that the motor speed exceeds 
the synchronous speed momentarily and settles down to the 
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Fig 441 Steady-state characteristics of a voltage 
coot rolled 3 phase induction motor of 
slip = 0,0311 and alpha = 80° 
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steady state sub-synchronous speed after some oscillations 
The start-up transient current with a load torque of 
3 N-m is shown m Pig 4 9 for a triggering angle of 70° 

Computed results of phase current i 1a> phase voltage v as and 
electromagnetic torque T 0 during steady state for a chosen 
triggering angle and load torque are shown in Pig 4 10 
These waveforms have very close agreement with the results 
given m Pig 2 8 of Chapter 2 It will he interesting 
to observe from Pig 4 11 that the rotor flux linkages (* 2 a 
and 4 > 2 q) a:ce practically constant during steady state opera- 
tion although the stator current contains large harmonic 
components The sixth harmonic components m the electro- 
magnetic torque and the transformed stator currents (i-j d 

and i 1 ) increases with increase in hold off angle Y • The 

* Q. 

details of the test machine are given in Table 2,1 of Chapter 2, 
4 3 2 Example 2 Digital Simulation of a Delta Connected 
Induction Motor 

The schematic diagram of a delta connected induction 
motor with voltage control is shown in Pig 4 12 The machine 
windings are permanently closed and a pair of back-to-back 
connected thyristors are included in each of the three lines 
The different possible modes of operation for such a scheme 
during transient and steady state are same as m the previous 
example and listed m Section 4 5 11 
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It may be observed that all the phase •windings of the 
motor would be carrying currents irrespective of the mode of 
operation This is because of the delta connection of the 
motor winding* During mode-5 operation the current flow is 
due to the presence of possible third harmonic voltages 
Thus it is obvious that this example problem is different 
fiom the previous example in respect of analysis and per- 
formance 


4 3 2 1 Formulation 

The equivalent circuit of a delta connected induction 
motor as obtained from Fig, 4 2 is shown in Fig 4 13 in 
which the possible mesh currents are marked During mode-1 
operation, there exist three independent mesh currents while 
during modes 2,3, and 4 operations only two independent mesh 
currents exist There is only one mesh current (i Q ) during 
mode-5 operation The general expression relating the wind- 
ing voltages and currents (refer equation (4*12)) is 


L" * 


. T 

~s 


p i s + R i s - I - L " p I 
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(4.54) 
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1b Ho* 


l*1a I 1'b I 1c * 


The objective is to obtain state equations corresponding to 
each mode of operation in terms of the mesh currents and source 
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voltages This can he easily achieved utilising the mesh 
analysis given hy Kron 15 3] 

Mode-1 

The mesh currents considered during this mode of opera- 
tion are i Q , i-| and ig The relationship between the phase 
currents and the mesh currents are 


-old “ G 1 -new 


(4 55) 


where 


i T 

—old 




1 1b 1 1e * * 

1 0 

0 1 

0 0„: 


—new " ^ 1 o 


X 1 


1 2l 


The matrix G 1 can be written by inspection from Pig 4 13 
The state equations (m terms of the new currents) are 


°1 K °1 P inew * -°1 E °1±new - °1 L " P I + °1 - 


(4.56) 


where 


T 

C 1 I 


f a'n " v b'n 
Vn “ v o'n 


Mode- 2 

Line a'a is open and hence independent mesh currents 
during this mode of operation are i Q &nd- 12 
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The connection matrix C 2 between old and new currents are 

1 0 

i 

(4 57) 


1 1 
1 0 


The state equation is 
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(4 58) 
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Mode-3 


In this mode of operation, the line b'b is open circuited 
The mesh currents are i o and 1 ^ where path of 1 ^ is indicated 
by dotted line (Fig 4 13) The connection matrix C, is 

r~ — f j 

1 0 


c 3 = 


1 0 

1 -1 


(4 59) 


The state equation is 
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Mode-4 


Line o’ o is open circuited The independent mesh currents 
are i and i-j The connection matrix het-ween l^ a » i>I 0 and 


the mesh currents are 

~1 1 

°4 = 1 0 


(4 61 ) 


The state equation is 


°4 C 4 P ine W * ' 0 4 R °4 ±new ' °4 ^ °4 1 (4 62) 


where 


li„ ii I » 


1 = [1 1j 

—new 1 o 1 




v a'n “ v h’n 


Mode-5 


Only one mesh current (i Q ) exists during this mode of 


operation The connection vector _c^ is 


and the state equation is 


°5 L s °5 p ^new 


-c T R c c l - Cr 1" p I c? 7 
n —5 new —5 — — 5 ~ 


(4 65) 


(4.64) 


where 
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4 3 2 ? Digital simulation, and results 

The formulation described m Section 4 3 2 1 is used 
for the digital simulation. Transient and steady state cha- 
racteristics are given m Pig 4 14 to 4 17 Starting torque 
as a function of time is shown in Pig 4 14 for ot * 90® and 
load torque T^ = 0 35 U~m Figure 4 15 shows the speed versus 
time curve The motor settles down to normal speed after a 
few oscillations The steady state characteristics are shown 
m Pig 4 16 and 4 17 The stator phase current is continuous 
as expected The electromagnetic torque (Pig 4 16d) has a 
predominant sixth harmonic component The rotor flux linkages 
if» 2 d and 4»2q are P rac-tlcall y constant (Pig 4 17) although the 
stator currents i^ and have high harmonic components 

433 Discussion 

Por chosen values of a and T , digital simulation re- 
sults of examples 1 and 2 are obtained for steady state ope- 
ration The slip, ♦ and average values of ^2q’ 1 1d 

and i^ a over one cycle are tabulated (Table 4 1) In order to 
facilitate a comparison between these results, same voltage 
per phase is taken for the wye and delta connected motors 
The results, given m Table 4 1, show that the characteristics 
of these drives are not identical especially at larger values 
of a The difference in characteristics m both cases is 
probably due to the presence of circulating current m the 
delta connected stator during single phasing Por smaller 
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Fig 4 14 Starting Torque 



Fig.447 Steady stale computed results for a 
A connected motor at ot=7Cf 3*04057 
(q) d-axis current 
(W q-axis current 

(c) d-axis rotor’ 

(d) : q-axis rotor 







1 50 6 0 0 0260 -1 087 -0 077 -2,922 2 757 2 311 44 48 

2 50 8 0 0.0414 -0 985 -0 185 -3 229 3 749 2 856 33 4 

■y 50 12 0 0 0969 -0 722 -0.349 -5 242 6 369 4 760 31 14 

/p 60 4 0 0 0164 -1 123 -0 020 -2 865 1 326 1 960 55 16 

5 60 6 0 0 0320 -0 968 -0 185 -2 935 2 870 2 370' 34 86 

6 60 8 0 0 0550 -0 814 -0 303 -3 693 3 892 3 090 29 93 

7 70 6.0 0 0471 -0 751 -0 313 -3 315 2 873 2 530 31 16 

8 70 8 0 0 1065 -0 487 -0 391 -5 410 4 402 4 030 34 62 

9 80 400 0358 -0 705 -0,290 -2 717 1 878 1 910 33 99 

10 so 6 0 0 1740 -0 258 -0 334 -6 500 3 859 4 360 41 73 

11 90 2 0 0 0153 -0 798 -0 163 -2 080 1 024 1.340 42 26 

12 90 3 0 0 0534 -0 420 -0 335 -2 571 1 547 1 730 37.88 

13 1 00 1 0 0 0068 -0 878 -0 3 27 - 2 086 0 6 86 1 260 5 4 1 7 

14 100 2 0 0 5897 -0 028 -0 099 -5 744 1 788 3 470 53 80 


contd . * 

f 

( 


f 



TABLE 4 1 (continued) 


i 

I 

1 10 


SI 

Eo 

a 

Load 

torque 

(W-m) 

1 

50 

6 

0 

2 

50 

8 

0 

3 

50 

12 

0 

4 

60 

4, 

,0 

5 

60 

6 

0 

6 

60 

8 

0 

7 

70 

6 

0 

8 

70 

8 

0 

9 

80 

4 

0 

10 

80 

6 

0 

11 

90 

2 

0 

12 

90 

3 

0 

13 

100 

1 

0 

14 

100 

2 

0 


a ^ or winding 

sup “ 

0 0260 -0 981 0 478 -1 169 3 830 2 310 44 74 

0 0413 -0 946 0 332 -0 91 2 4 870 2 851 33 46 

0 0965 -0 802 0 059 -1 340 8 110 4 722 31 05 

0 0163 -0 983 0 544 -1 570 3 02o 1 962 55 30 

0 0331 -0 936 0 314 -1 051 3 906 2 341 33 51 

0 0550 -0 858 0 144 -1 257 5 197 3 082 29 32 

0 0470 -0 808 0 103 -1 424 4 133 2 520 31 04 

0 1060 -0 620 -0 096 -2 463 6 436 4.011 34 50 

0 0359 -0 757 0 098 -1 400 2 969 1 889 33 80 

0 1687 -0 399 -0 161 -3 633 6 490 4 290 41 40 

0 0172 -0 776 0 236 -1 170 1 306 1 244 40 30 

0 0544 -0 529 -0 088 -1 430 2 607 1 738 37 68 

0 0067 -0 779 0 411 -1 463 1 598 1 253 54 10 

0 5436 -0 080 -0 077 -3 926 4 357 3 383 53,48 


Ill 


values of Y the difference m characteristics is not appre- 
ciable The differences m the d and q axes current compo- 
nents between the two cases are due to the different references 
adopted m calculating these quantities Hence, only the r m s 
current magnitudes can be compared 

Although the problem of the analysis of a voltage con- 
trolled delta connected induction motor has been mentioned 
in reference 122] , a comprehensive analysis and simulation 
has not been attempted in the literature so far 

4 4 CONCLUSION 

Two models of an induction motor are described m which 
the first is particularly adapted for the simulation of motors 
where stator terminal switching is employed for speed control 
application The second model is well suited for rotor ter- 
minal switching conditions The ma^or features of these 
models can be summarized as follows 

1) The time dependence of the inductance matrix is 
eliminated 

2) Identity of the motor terminals, where switching is 
employed, is retained 

3) Simple but exact equivalent circuits are provided for 
either stator or rotor windings so that foimulation 
of the problem with any arbitrary external network 
is easily obtained with well-known methods of network 


analysis 
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The major contribution in this chapter is the develop- 
ment of equivalent circuits for stator and rotor windings 
of the induction motor which are valid for both dynamic and 
steady state conditions 

The transformation of variables anto synchronously 
rotating reference frame has certain important significance, 
not only for sinusoidally excited machines, but also for 
thyristor controlled drives Simulation results indicate 
that the rotor flux linkages are essentially constant even 
•when the stator current contains large harmonics due to 
thyristor switching The constant nature of flux linkages 
leads to two important features of the model presented m 
this chapter 

1) Tnc constant t>ait of the deieloned toraue (eQn 4 1 q ) 

(the torque contains a constant term with harmonic com- 
ponents m all thyristor controlled machines) is pro- 
duced by the fundamental component oi stater current 
alone This property i p aoulicnble not only to voltage 
controlled machines but -tUo to many other tourist or 
controlled drives such as volt age-fed m/erters 
2) The constancy of rotor flux linkages leads to certain 
simplified analysis whicn nemits the development of 
transfer function between tne torque and the control 
variable such as a This is particularly important 
m the design of sneed controllers This aspect is 
considered m the next chapter 
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The f omul at ion and simulation for two different 
conf igurations of voltage controlled induction motor are 
given to illustrate the applicability of the model Computed 
results of transient and steady state characteristics are fur- 
nished for each example In particular, the example of a vol- 
tage controlled delta connected induction motor has not been 
attempted m the past and the results of this example are 
compared with the results obtained for the wye-connect ed 
mot or 


CHAPTER 5 


SIMPLIFIED ANALYSIS OP THE STEADY STATE PERFORMANCE 
OF A VOLTAGE CONTROLLED INDUCTION MOTOR DRIVE 

5.1 INTRODUCTION 

The previous chapters considered steady state performance 
of a voltage controlled induction motor using digital simula- 
tion While this may he adequate for the analysis of a motor 
with open loop speed control, it is not sufficient for the 
analysis of a speed control system with feedback Vhat is 
essential is a method ^o directly relate tne motor performance 
to the control variables, (sn~h a*- firing angle, source ■'oltage) 
Thi^ chanter credent a simplified analysis based on the hybrid 
model t_ v en in Cnapter 4, tj calculate t^e steady sta + e per- 
formance of a voltage controlled induction motor hhile 
otner conj igurat ions can also be considered, the analysis 
given here is apolicable to a wye-connect ed motor with a pair 
of back-to-back thyristors connected m each line 

It was observed from the results of simulation given m 
Chapter 4 that the rotor flux linkages referred to a synchro- 
nously rotating reference frame are essentially constant even- 
though the stator currents contain appreciable amount of har- 
monics This observation can also be confirmed analytically 
and is done later in this chapter This is an important 
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observation because it justifies the useful approximation, 
of neglecting harmonics in the rotor flux linkages With this 
approximation and examination of the torque eauation given 
m Chapter 4, it can be shown that only the fundamental com- 
ponents of the armature currents give rise to the constant 
component of the torque developed This fact is utilized 
m this chapter to derive three nonlinear algebraic equa- 
tions from the solution of which, one can calculate the 
steady state performance of the motor for a given slip* 
firing angle and the source voltage. 

The hybrid model given m Chapter 4 1° used to f lr'-t 
derive tne expre Sion’ 7 * for the stator current^ and terminal 
voltage of tno motor The calculation of the current leads to 
a nonlinear algebraic eauation m teims of + he p^ase angle (4*) 
ana the rotor flux linkages (4*^ and ^q) n he r of° r flux 
linkages are solved m term 1 - of the fundaneru- 3 ! conronents of 
Ire stator currents *hicb m turn are calculated from tne 
fuid mental components of the terminal /oltagc T^is leads to 
s set of two additional algebraic equations ~Y c combined set 
oi nonlinear algebraic equations ere solved using A <ewton- 
Raphson method to calculate the steady 'fate performance 
of the motor The slip (5), the tnggcnrg sr rle (a), and 
source voltage magnitude (V^) appear m the e aunt ions ns para- 
meters The fundamental current and voltage magnitudes, phase 
angle ($) and the torque can be calculated in terms of a ,5 and 
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Two examples are presented to illustrate the analysis 
and its application 

5 2 SIMPLIFIED STEADY STATE ANALYSIS 

521 Assumptions 

The major assumption used here is that the rotor flux 
linkages re f erre & to a synchronously revolving 

frame are constants during steady state operation from the 
simulation results given in Chapter 4 (Fig 4 11), it can he 
observed that these flux linkages are constants even if the 
armature currents are nonsmusoidal due to controlled swit- 
ching The validity of this Qr,c, u^rtion is al°o "proved snn- 
lyuically m a later section (5 4) Rofeinnf to the torque 

cx wcs^ion (T * n (i 1d ^ 2 q ” X 1 o li can be seen 

th^t tne coi°tont component of toroue is produced by the d c 
coiij onents of i 1d and i<j q since ^d anii ^2q n ' re constants The 
fund mental component of stator current contribute** the d c 
conponenis of ana i«j ^ 

The second assumption is that t c analysis is restricted 
to inode-1 operation (y < 60°) in wnich either t^ree or two 

■windings of the motor are connected to the pow-n? source Tne 
st-’tor windings are assumed to he wye connected /ith isolated 
neutral 

5 2 2 Analytical Expression for the Stator Current 

Typical voltage and current waveforms during mode-1 
operation are given m Fig 5 1 ^or convenience, the half 
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cycle period is divided into six sections starting from 
w t = a as shown m Fig 5 1 Referring to the circuit model 
(Fig 4*3) the differential equation defining the stator 
current i^ a can be written as 

R 1 1 1 a + L 1 p( ' 1 1a + I 1a^ = V a’n ” v s = Y as ^ 1 ^ 

where v g is the voltage at the star point of the motor with 
respect to the source neutral, v , is the source voltage and 

cl il 

v is the voltage across the a-phase winding From the cir- 
as 

cult model (Fig 4 3) v_ can be calculated .using Millman's 

s 

Theorem 154], for different modes of operation, as 

v g sa 0 when all phases conducting (5 2) 

v g = ~ p l^)/2 when phase b non-conducting (5 3) 

v o = -(v c , n - L}' p I 1f> )/2 when phase c non-conducting (5 4) 

Since and ^ 2 q are constants during steady state, the 
time derivatives of the dependent current sources are given by 
(see equation (4 16) of Chapter 4) 

P I 1a - E (- *2d Sin 6 o + *2q cos 0 o )/L 1 
p I 1b = K{- sm(e o -2V3)+<f> 2 q oos ( 6 0 ~ 2 V 3 ) >/II^| , 

p I 1c - K { - tf> 2a sin(e o +2ir/3)^ 2q cos (0 o +2V3)tyl»!j 

where 


(5 5) 
(5 6 ) 
(5 7) 


K 


« l 12/ L 2 
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The solution, of equation (5 1) for each section can be 
easily obtained as it is equivalent to a R-L circuit excited 
by a sinusoidal voltage with known initial condition The 
current expressions as derived in Appendix F are 

section 1 

i i a = pnjl {sin (0 o + + °l ex P(~ e(e 0 - a)) } (5 8) 

section 2 

i 1a 55 t§YT f ^ 3/ 2 sm(e o+ ir/6+5) + C 2 exp(-e(e o -*/3-4> )) (5 9) 

section 3 

1 1a = * sin (9 0 +^) + °3 exp(- $ (e o -*/3-o )) } (5 10) 

section 4 

1 1a = fin 1/3/2 smto 0 -*/6+«) + c 4 exp(-B(fl 0 -2*/3-*))) (5 11) 
section 5 

i 1a = lifl * sin ( 0 o + ^ + °5 ex P(~P<V 2 * /3 “ «))) (5 12) 

section 6 

i i = 0 (5 13) 

1 a 

where 

E ={ E^+Eq>^, E a * -K i|> 2q » E q 85 V m + K ^2d» Z 1 = R 1 + 3 X 1» 
X!(*a>L!j, 0 * Hj/XJ, P - tan^O/*), «- 6 - *> > 

0 = tan" 1 (%/ E q)» 
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C-j = a c 6) 

C 2 = C 1 exp(- $( ir/3+«J-a)+siu(ir/3+<jH-6)-/5/2- ros(* + 6) 

V 

{ 

Oj = 0 2 exp(~ B(a-<J>)) + /3/2 cos(a + 6) - sm(Tt/3+a+fi) 

, = 0^ exp(-B(ir/3+<}>-a) )+sm(2ir/3+<f»+6)-» / 3/2 cos(<{>+6) 

and 

= C 4 exp(- )+^3/2 cos(ct+6)_ sm(2ir/3+a+6) 

At the end of section 5 (0 = tt +<t) , 1 - 1 « 1360011163 zero 

o * £L 

Equating i^ a to zero at this instant, we get a nonlinear 
equation defining the phase angle <{> as 

0-j.j exp( -Ba) + C^ 2 exp(- B<J>) =0 (5 14) 

where 

= - / 3/2 cos(<t> + 6) {1+exp(-0ir/3) } -sm(4> +5) exp (Bir/3) + 
sin(27r/3+<{>+6 ) + sin(ff/3+<J>+ 6) exp(-B^/3) 

C^ 2 =*^~3/2 cos"(a+6){ 1+exp(- B* /3) ) -sin( a +6) exp(-B2ir/3) - 
sin(2ff/ 3+CH-6) - sm(ff/3+a+fi) exp(-Btr/3) 

The solution of equation (5 14) for an operating point (a, 

S) cannot "be obtained as the term 6 is a function of 'l>2d and 
$ 2 q A l inear relationship between t|> 2( i> ^ 2q x ia* X 1 q 
during steady state is obtained from equation (4 17) by setting 
p ij> 2d and p i|> 2q to be e< 3 ual to zero The relationship is 

given "by 

% 

* 

t 


t 


f 
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1 1d 

1 

n 2 s x. 


*24 


Ro m 




_ 1 1 q_ 

C 

-s x 2 r 2 


2q 


where 

m as ^3/2 Xg ss oj 

5 2 3 Analytical Expressions for the Stator Voltage 


(5 15) 


sectional, 3 and 5 (Pig 5 1) all the three phases 

conduct and hence, v a m equation (5 1) is zero Therefore, 

s 


v = v , 
as 


. = V sin 6 

a'n m o 


(5 16) 


During Section 2, c-phase is non-conducting while a and h 
phases carry current The a-phase voltage is given hy 


Replacing, v g hy equation (5 4) and eliminating the derivative 
of I 1o hy equation (5 7), we get 

v = H ^3 T sin(e +*/6)+Bi sm(e +2ir/3)*E ? cos (e +2^/3)) 
as ^ m 0 1 0 ^ o 

(5 17) 

where 

E 1 - K ♦ 2a , E 2 “ K 41 2q 

During section 4, only a and c phases conduct while h does not 
conduct The neutral voltage v g is defined in equation (5 3) 
during this mode The stator voltage is given "by 


122 


~ i (*^3 sin(0 o ~Tr/6)fE^ sm(e 0 -2 tt/ 3)*E2 cos(8 0 ~2tr/3) } 

(5.18) 

In section 6, a-phase is open -circuited A back-emf is 

induced in the a-phase winding due to presence of rotor currents 
and is given by 

v = Li 1 p L 
as 1 r n a 

» -E.j sin e Q + E 2 cos e Q (5 19) 

5 2 4 Eouner Analysis of the Voltage 'Waveform 

Since the voltage waveform is completely defined 
(5 16 - 5 19) » it is possible to derive analytical expressions 
for the Fourier components The even harmonics and d c 
componaits are absent m the voltage waveform due to half 
wave symmetry The fundamental Fourier components of stator 
voltage as derived in Appendix G- are 

a 1 = Tfz t ( V m +E 1 ^ ( G0S 2ot ” cos 2 ^ + 

E 2 (sm 2 a - sin 2<|>+2a-24>)} m /z/3 v^ (5.20) 

b 1 * { -E 2 (cos 2 a - cos 2<j> ) + V m + 

( V Ei)(sm 2a -sin 2<f> -2a + 2$) « /2/3 v go ) (5 21) 

Expression for and b^ are given m Appendix G 

In terms of the fundamental components of voltage and 
current, the stator side equivalent circuit can be represented 
on a per phase basis as given in Fig 5 2 Referring to Fig 5 2, 
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an expiession relating the rotor flux linkages to the current 
and voltage can he derived as follows 

let , 

1 = x do “ 3 i qo ’ Y = V do " 3 V q 0 » X = X do “ 3 X qo 
Referring to equation (4 16) of Chapter 4, the dependent current 

A 

I can he related to the rotor flux linkages as 

5 - 'V* ipf- <*2ao - 3 <5 22) 

From the circuit model (Pig 5 2), the expression for i is 
given by 

, v - d I X 2 

i = am (5 23) 

"1 


where 

X!j = w Ii!j i Z!j = R 1 + 3 Xj 


After simplification, we get 


^o ' 3 X q 0 * 


< Cv d - /V2 K 4> 2d )-3(v„ +•'1/2 K + 2q )} (R-,-3 Zp 


l«? 


qo 

T 


(5.24) 


equating real and imaginary parts of equation (5 24) and 
expressing in the matrix form, we get 





+ 


*3/2 K 


-h -a 


*240 

a -h 


'‘‘290 


(5 25) 
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where 

a = R-,/1 Zlj'j 2 , b = V\/\Z'\\ 2 

It may be noted that the induction motor model shown in 
Rig 5 2 is equivalent to the conventional passive equivalent 
circuit for sinusoidal excitations m steady state This can 
be proved (see Appendix H) by showing that equation (5 23 ) is 
equivalent to the equation defining passive equivalent circuit 

5 2 5 Prediction of Operating Vector 

The operating vector of a voltage controlled induction 
motor consists of phase angle <f> , fundamental components of 
current and voltage and constant component of torque The 
operating point is defined by triggering angle a , slip S and 
the source voltage magnitude V m The problem is to predict 
the operating vector for a given operating point, An analyti- 
cal expression for <f> is given in equation (5 14) During 
steady state, ^2q are ^nearly related to 1 ^ and i^ q 

by equation (5 15) Since the harmonic components i-j^ and 1 ^ 
have very little effect on and ^q* it is sufficient to 
know the fundamental components of stator current to determine 
^2d ^2q & urin 6 steady state The fundamental component of 

current results from the fundamental component of applied vol- 
tage The current and voltage are related by the parameters of 
the passive equivalent oircuit as 




» 

ii 


(i do " 3 V (R eq + 3 X e<t> 


(5 26) 
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•where R eq and X eq are the equivalent resistance and reactance 


of the pass^y^ pqurvalent circuit, given by 

X? 

R . * L + ? S _£ 

69 ( V s > + - * n > 


(5 27) 


and 


‘eq 


( ya) ♦ + V 

Sl (R 2 /3) 2 7 + X rl ) 2 


X 


m 


(5 28 ) 


•where R^ , R 2 » x s i» x r p and ^ are -fclle parameters of conventional 
passive equivalent circuit 


On separating the real and imaginary parts of the equation 


(5 26), we get 










V, « 
do 

B 1 do 

R eq + 

1 qo 

X oq 





(5 

29) 

V 

— I 

R 

la 

X „ 





(5 

30) 

qo 

qo 

eq 

do 

eq 







Combining 

equations 

(5,29) and 

(5 30) 

with 

(5 20) and (5 

21) 

> 

we get 











r? 4tt 

v 5 y 

^do 

R on + 
eq 

1 qo 

z eq> - 

(V + 
in 

K *2d 

)(cos 2a 

- cos 

2<t> 

) - 


K i|» 2q 

(sin 2 

a - 

sin 2<J> 

+ 2 ot - 

2* ) 

= 0 


(5 

31) 

72* 4 it 

h t 

(l qo 

R eq - 

x do 

V + 

K *2q 

(cos 

2a - cos 

24> ) 

- 



(V +K ^ ) (sin 2 a*. sin 2<j> - 2a + 2 4> ) - y v m = 0 

(5 32 ) 

Making use of (5 15) » ‘the nonlinear equations (5 14) » (5 31) 
and (5 32) w <J> , x dQ and i qo are solved simultaneously by 


126 


Newt on -Rap hs on. method { 55 ] for a given operating point It 
may be remembered that the equations derived in Section 5 2 
are applicable only for mode 1 operation ( y < 60 °) Hence, 
the mode of operation is to be checked during the solution of 
the nonlinear equations A flow diagram giving the details 
of mode identification and the solution is shown m Fig 5 3 
If the hold off angle y is greater than 60 °, the programme 
returns without computation Proper initial values ( l^ 
and iq x ) which are generated using passive equivalent circuit 
parameters, are provided for faster convergence The Jacobian 
which consists of the derivation of the three functions with 
respect to <J> , i^ 0 and ig 0 are easily calculated because of the 
availability of analytical expressions for these functions 
Once the steady state solution is obtained the torque is 
calculalcd using equation (4 19) The expressions for the 
r m s values of current and voltage are given by 


1 t 2 2 

*do + 1 qo J 

(5 

33) 

1 r 2 2 1 * 

-tv, + V _ * 
do qo 

(5 

34) 


5 3 EXAMPLE 
5 3 1 Computed Results 

The steady state performance characteristics of the test 
machine (Table 2 1) were predicted using the simplified analysis 
discussed m the earlier Sections Predicted values of $ , 
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Fig 5 3 Flow diagram ^predict operating vector 
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*2d* *2g> lyjjjg and fundamental power factor (cosine 

of the angle between fundamental components of voltage and 
current) for mode-1 operation are plotted against slip S m 
Pigs 5 4 to 5*7 It can he observed that all these varia- 
bles change rapidly for small changes in slip when the motor 
is near maximum speed This Is because of the fact that the 
phase angle $ changes rapidly (Pig 5 4) near the maximum 
speed for small variations m slip resulting in large varia- 
tions m the hold off angle y Por larger values of slip 
(S y 0 2), the variations are not very significant Varia- 
tions of the rotor flux linkages for different values of a 
are plotted against slip S in Pig 5 5 The magnitude of 
the fundamental components of voltage and current are given 
m Pigs 5 6 and 5 7 In Table 5 1, predicted and simulation 
results (4> , and ty 2q ) are furnished for comparison 

In order to check the accuracy of the predicted results, 
the steady state characteristics of a motor described in. 
reference 114] are computed The computed torque-speed 
characteristics of this motor is plotted in Pig 5 8 for 
constant values of ct for comparing with the published re- 
sults { 14J The parameters of the motor are given in 
Table 2 2 
5 3*2 Discussion 

On comparing with the results published in [141 , it 
can be seen that the predicted torque-slip characteristics 
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COMPARISON BETWEEN PREDICTED AND SIMULATION RESULTS 


CO 

p 

H 


0 

PI 

a 

o 

H 

-e 

0 

H 

§ 

H 

tn 


CO 

p 

H 

2 

to 

0) 

PS 

rd 

CD 

P 

o 

H 

nd 

CD 

P 

P 


CTM 


(No 

"3- S 


CO 

CD 

CD 

u 

W) 

CD 


wjP 

O 1 | 
cm A 


'd+j* 


CO 

CD 

s 

(D 

nd 


P _ 
H pS 

rH 

CO P 


CO 

03 

CD 

U 

w 

0 

fd 


H O 
CO 


Cn 

£> 

O 

in 

CO 

VO 

CO 

CM 

CM 

tn 

VO 

CM 

o 

CM 

VO 

T*** 

CM 

T~* 

cn 

CO 

tn 

CM 

cn 

CM 

cn 

T““* 

n- 

cn 

O 

-d- 

tn 

o 

1 

o 

1 

O 

I 

o 

i 

o 

1 

o 

1 

O 

I 

o 

i 

o 

t 

o 

1 

t 

o 

l 

o 

i 

VO 

0D 

o 

tn 

co 

cn 

o- 

in 

00 

CM 

CM 

Cn 

T- 

in 

CO 

cn 

cn- 

cn 

<• 

in 

in 

CO 

cn 

o 

*c~ 

018 

o 

cn 

VO 

v~ 

i 

o 

i 

o 

l 

O 

* 

1 

o 

1 

o 

i 

O 

i 

o 

l 

O 

t 

V* 

1 

o 

i 

o 

C\) 

oo 

tn 

CO 

00 

o 

CM 

o 

t>* 

o 

00 

T"* 

V“ 

T~ 

CM 

o 

cn 

CM 

■sf- 

^d" 

cn 

m 

cn 

CM 

o 

«d- 

in 

in 

cn 

cn 

CM 

tn 

cn 

«d* 

tn 

in 

cn 

*d" 

V“* 

cn 

CT\ 

£> 

O 

o 

00 

V* 

00 

00 

CM 

tn 

cn 

t- 

CM 

O 

in 

vo 

r~ 

vo 

T* 

m 

00 

cn 

cn 

CM 

cn 

cn 

cn 

O 

00 

* c d" 

cn 

O 

1 

o 

1 

O 

1 

o 

i 

o 

1 

O 

1 

o 

i 

O 

i 

O 

1 

O 

o 

\ 

o 

i 

086 

cn 

00 

cn 

o 

tn 

CO 

tn 

CM 

■<d- 

cn 

T*" 

in 

CO 

cn 

00 

cn 

o- 

•<d- 
<• 
l n 

186 

*d" 

O 

00 

o 

o 

D— 

vo 

V“ 

1 

o 

i 

o 

1 

O 

\ 

V 

1 

o 

i 

o 

1 

O 

1 

o 

1 

O 

1 

T~ 

1 

o 

l 

m 

r- 

*T“ 

00 

VO 

in 

cn 

CM 

cn 

vo 

in 

cn 

CM 

oo 

00 

T*“ 

o 

in 

cn 


cn 

cn 

cn 

CM 

cn 

cn 

LfS 

in 

cn 

cn 

CM 

cn 

cn 

cn 

<■ 

rn 

in 

m 

o 

cn 

VO 

CM 

O 

T** 

O 

■*d" 

vo 

o 

cn 

cn 

CM 

t> 

T" 

o 

V* 

cn 

o 

00 

in 

o 

«a- 

CM 

t— 

' t d~ 

in 

cn 

vo 

cn 

vo 

CM 

O 

vo 

in 

o 

o 

O 

o 

O 

o 

o 

o 

O 

o 

o 

O 

o 

o 

LTN 

o 

in 

o 

in 

o 

in 

60 

60 

60 

09 

60 

60 

O 

tn- 

o 


T -<MtO»^^ VDt ^' 00 Cr '®P^ 


134 


( 


»j ' 

I jt 

I 

«* f 





If 

u 


IF _ 

gpr 




I co 

I a) 

p $4 
I tuO 
I 0) 

I d 


o 

o 

in 

VO 


00 

CM 

tn 



tn 

CO 

VO 


t> 

t> 

in 

tn 

in 

CM 

tn 

t> 

c n 


O 

^r 

£>~ 

tn 

tn 

CM 

o 

V* 

tn 

tn 

CM 

T~ 

o 

tn 

tn 

CM 

x — 

T— 

o 

i 

o 

o 

O 

o 

O 

O 

o 

o 

o 

o 

O 

o 

o 

l 


l 

i 

l 

1 

1 

1 

i 

i 

1 

t 

1 

CM 

o 

in 

CM 

CM 

00 


o 

CO 

CO 

T~ 

in 

cr\ 


C^ 

o 

VO 

CM 

VO 


tn 

cn 

00 


o 

tn 

tr- 

in 

tn 

CM 

-r- 

Ov 

in 

CM 

V“ 

o 

00 

tn 

CM 

r- 

O 

O 

o 

o 

T- 

o 

O 

O 

O 

o 

o 

o 

O 

O 

O 

O 

t 

1 

1 

1 

i 

l 

1 

i 

1 

\ 

\ 

1 

1 

t 

CT\ 

00 

V 

in 

cr\ 

tn 


in 

Lf\ 

<r 

r~ 

tn 

tn 

c0 



T- 

V 

CO 

CM 


in 

o 

CO 

tn 

o- 

C' w * 

C— 

CO 

VO 

o- 

T~ 

tn 

CM 

00 

in 

VO 

CO 

tn 

VO 

o 

CM 

tn 

xtfc 

o- 

x}- 

tn 



un 


tn 


*M“ 

in 

in 

tn 

tn 

in 

in 

T~ 

cr\ 

tn 

vo 

c 

c— 

in 


00 

tr~ 

r>~ 

tr*- 

in 

tn 

VO 

CM 


t— 

cn 


O 

m 

c- 

tn 

tn 

CM 

o 

V* 

tn 

tn 

CM 

v* 

o 

tn 

tn 

CM 

T ~ 


o 

t 

O 

I 

o 

o 

i 

o 

i 

o 

i 

o 

1 

O 

o 

1 

o 

o 

i 

O 

1 

o 

1 

o 

i 

!P 

i ^ 

o 

o 

CM 

in 

VO 

T“~ 

■*4“ 

CM 

CTn 

cr\ 

in 

in 

o 

CM 

in 

tn 

cn 

o 

o 

Ov 

in 

tn 

cn 

cn 

V* 

tn 

tn 

CTs 

tr- 

Ct 

in 

o 

19 

O 

1 

r~ 

1 

O 

l 

o 

i 

O 

1 

o 

1 

o 

1 

o 

\ 

O 

i 

O 

1 

O 

l 

O 

{ 

o 

1 


fm ; 

1 

1 

CO 

CD | O 

CD | ^ 

o 

tn 

O 

20 

o 

c- 

cn 

v- 

VO 

tn 

(T> 

O 

tn 

t- 

00 

CM 


CM 

t> 

1 U 1 

t I ou 

I cu i tn 

vo 

c- 

o 

-r* 

tn 

tn 

CM 

CO 

V~ 

in 

C^ 

-M- 

0) 

tn 

tn 

v£> 

O 

in 

i 

i 

d | 

I 

l 













✓—v 

p* . 
H 2 

i 

1 

1 

t 

i 9 

1 VD 

i t™ 

O 

in 

tn 

o 

o 

c n 

V" 

O 

in 

in 

O 

CO 

T*“ 

CM 

f“ 

tn 

O 

c- 

in 

tn 

r- 

O 

T~* 

tr- 

o 

in 

-Mr 

*r 

CM 

tn 

CM 

vo 

T~ 

! * _ 
CO 

\ O 

1 

o 

o 

O 

O 

o 

o 

o 

o 

o 

O 

O 

O 

CQ 

0) 

1 

1 

1 













o 

<D 

rd 

\ 

1 

1 o 

1 r- 

1 

70 

80 

80 

80 

80 

80 

08 

06 

06 

06 

06 

06 


ITN 


QO CTi O 


CM tn 


in vo 


& Cvl W CXJ CM CVJ CM 


136 


(Fig 5 8) have very good correlation with the exact values 
This validates the correctness of the assumptions made in 
the simplified analysis According to the induction motor 
model developed in Chapter A (equation 4 19), the constant 
component of the torque is produced by the fundamental com- 
ponent of stator current primarily because the rotor flux 
linkages referred to the synchronously involving axes are 
practically constant 

The predicted values of , 4> 2d , for different 
values of oc and S are shown m Table 5 1 for the test 
machine Results obtained from digital simulation (Chapter 4, 
example 1) are also given m the same table to enable a com- 
parison between the two results Very good agreement between 
the two results is evident from Table 5 1 

5 4 nFFECT OF STATOR CURRENT HARMONICS ON ROTOR FLUX 

LINKAGES ij> 2d AND # 2 

It is stated in Section 5 2 that and $ 2 q 3X6 

essentially constant during steady state and harmonic com- 
ponents of stator current have very little effect on them 
The digital simulation results illustrate the validity of 
this assumption An analytical expression to justify this 
assumption is derived m this section 

Consider the differential equations defining the rotor 
flux linkages derived in Section 4 2 1 of Chapter 4 given by 


p i[» 2d » ~ R 2 ^2d/ L 2 “ Sw ^2q + R 2 m Md^ 


137 


(5 35) 

P ^2q * " R 2 *2</ L 2 + Scu *2a + R 2 m i 1 q/ L 2 (5 56 } 

The harmonic components of ij> 2 d and m steady state result 
from harmonic components of i 1d and i 1q These can he calcu- 
lated as follows 

Multiplying equation (5 36) hy -j and adding to (5 35), 


pOJ>2d~3 ^2q^ * ~ R 2^2<r 3 ^2q)/^2" 3 S w ^ J 2d' 3 ^2q^ + 

R 2 m(i 1a -3 l 1q )/I 2 
In vector notation, 

A A A 

1 2 P ^ 3 -(R 2 + D S Xg) $2 + ^2 m 1 (5 37) 

where 

*2 = *24 -3 *2q> 


During steady state, operator 1 p ' can he replaced hy an w 
where n represents the order of harmonics and we get 


A 

(Rg+3 (n+S 5 X 2 ) ~ R 2 m ^n 


(5 38) 


where 

% 


2n “ "’Sdn " , 2qn- 


A 

x n - 1 


dn ~ 3 1 qn 


The expressions for i fln and 1 ^ in terms of the Fourier com- 
ponents of stator currents (see Appendix I) are 


! dn = ( a n-1 +a n+1 )2 + ( Vl + Vl )2} 

1 qn “ 2 i ^ a n-1“ a n+1 ^ + ^n-l^n+1^ } 


(5 39) 
(5 40) 


are 


where a n and b^ are the Fourier components of stator wrr-eat 
On rearranging the terms m (5 38), we get 
R-o m (Ro la — (n+S) £ 0 i„„) 


Rg + (n+S) T z 
Ro m(R o i„„ + (n+S) X, i, ) 


^ R z + (n+S) ir 2 

Equations (5*41) and (5 42) indicate that ^dn ^2qn 

are reducing rapidly as the order of the harmonics increases 
since the reactance increases from S X 2 for d c component 
(of i 1d and i 1q ) to (n+S) X 2 for the nth harmonic As an 
example, at a slip of 0 1 , the reactance for the sixth har- 
monic, is about 61 times greater than that for the d c com- 
ponent 

Using equations (3 41) and (5 42), * 2dn and «> 2an are 
calculated and the results are given m Table 5 2 along with 
the corresponding values of i dn and i qn It can be obser- 
ved from Table 5 2 that t 2dn ^ *2q» are a P^ roachin S zero 
for higher values of n Even for the sixth harmonic (the 
lowest order harmonic in i dn and i qn ) component of i 1d and 
i the magnitude of the flux linkages (harmonic) becomes 

quite small 

Fourier components of a n and b Q of the stator current 
can be obtained in three different ways The first method is 
to derive analytioal expressions for a n and b Q from the current 
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TABLE 5 2 

HARMONIC COMPONENTS OF E-Q CURRENTS AND ELUX LINKAGES 


a = 50° , slip =01 


n 

(& ) 

V“) 



0 

-5 390 

6, 410 

-0 711 

-0 350 

6 

0 845 

0 761 

-0 0012 

0 0014 

12 

0 233 

0 208 

-0 0002 

0 0002 

18 

0 025 

0 032 

-0 0000 

0 0000 


a =60°, slip = 0 1 

n 

i d (n) 

lq(n) 


V n) 

0 

-5 460 

5 370 

-0 6161 

-0 3831 

6 

1 101 

1 183 

-0 0019 

0 0018 

12 

0.090 

0 113 

-0 0001 

0 0001 

18 

0 102 

0 110 

-0 0001 

0 0001 


f 


1 


Table 5 2 (continued) 
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* = 70° , slip = 0*1 


n 

i a W 

i q {n) 

V»> 

V n) 

0 

-5 500 

4 330 

-0 4866 

-0 3939 

6 

1,186 

1 498 

-0 0024 

0 0019 

12 

0 119 

0 173 

-0 0001 

0 0001 

18 

0 085 

0 107 

-0 0001 

0 0000 




a = 80® , 

slip =01 

n 

i d (n) 

iq(tt) 

l| V n > 


0 

-4.855 

3 023 

-0 3867 

-0 3797 

6 

1 152 

1 5 42 

-0 0025 

0 0019 

12 

0 205 

0 357 

-0 0003 

0.0002 

18 

0 056 

0 078 

-0 0000 

0 0000 




a = 90® , 

slip = 01 

n 

i d (n) 

iq(n) 


*2q< n) 

0 

-4,115 

1 918 

-0 2684 

-0 3362 

6 

1 160 

1 406 

-0 0022 

0 0019 

12 

0.216 

0 399 

-0 0003 

0 0002 

18 

0.076 

0 179 

-0 0001 

0 0000 
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expressions (equations (5 8 - 5 13)) Another method is 
to make use of the Fourier components of voltages and harmonic 
equivalent circuits tl9J The third method is to obtain a Q and 
b^ numerically -while doing digital simulation The calcula- 
tions for the results given m Table 5 2 are performed using 
the first method and checked using the third method 

5 5 CONCLUSION 

The analysis presented m this chapter has two special 
features The first one is the faster computation of steady 
state characteristics for any given operating point Although 
the analysis is based on certain assumptions, the accuracy of 
the computed results is very good The computation time 
needed to predict the characteristics is approximately 1/50th 
of the time needed by state-variable method 

The second and most important feature of the analysis is 
that it can be used to study the dynamic performance under 
perturbations in load torque and supply voltage From the 
analytical expressions developed m this chapter, linearized 
model and transfer functions of the voltage controlled indu- 
ction motor can be developed which would help m the design 
and development of closed loop control schemes This aspect 
is considered in the next chapter 
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DYNAMIC ANALYSIS AND DESIGN OE A CLOSED-LOOP SPEED 

CONTROLLER 

6 1 INTRODUCTION 

In Chapter 5, analytical expressions for the steady state 
currents and voltages are derived and expressions for the phase 
angle and fundamental Eourier components of stator voltage 
and current are obtained as function of triggering angle, supply 
voltage and rotor flux linkages This enables the development 
of a linearized model of the drive and its associated transfer 
functions for the study of dynamic response under perturbations 
m load torque, triggering angle or supply voltage This can be 
used for the design of closed-loop control systems to regulate 
speed* 

This chapter presents the linearized state space model of 
the drive Stability of the drive, based on the eigenvalue 
analysis, is discussed Transfer functions of the drive are 
developed to predict dynamic response Using these transfer 
functions, a feedback control scheme is developed to regulate 
speed An example showing the design of the controller using 
classical control theory! 561 is given The response of the 
closed-loop system is predicted analytically and checked by 
experimental investigations 


6 2 LINEARIZED MODEL OE AN INDUCTION MOTOR WITH THYRISTOR 
VOLTAGE CONTROL 

In developing this model the stator transients are ignored 
and the motor is represented by a third order model given by 
(see Section 4 2 1 of Chapter 4) 


and 


^2d ~ ” R 2 ^2d/ L 2 ~ Sw ^2q + R 2 m 1 1d //l 2 

(6 1 ) 

P *2q " “ R 2 ^2q /L 2 + So » ^d + R 2 m 1 1 q/^2 

(6 2 ; 

p s - - (VV /(W J) + O-s) D / J 

(6 3; 

T q « m ^ 1 l d ^ 2q “ 1 1q^2d^ L 2 

(6 4: 


On linearizing equations (61-64) and arranging in 
matrix form, we get 




A\p _ 



2d 


2d 


1 1d 

A *2c 

= t 1 

**20 

4* B-4 



+ CAT, 


(6 5 
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The vector containing &i 1d and fii 1q in equation (6 5) 
is to he replaced in terms of the state variables AiJjg^, A<f>g q 
and AS The relationship between the fundamental components of 
voltage, current and flux linkages (see equation (5 13) of 
Chapter 5) is given by 



On linearizing equation (6 6), we get 



The fundamental components of the stator voltages (see Section 
5 2 4 of Chapter 5) are defined by 

7 a “ '! fi { <V K *2a )(coa 2a - cos 2$ ) + 

K t|» 2q (sin 2 a - sm 2<J> + 2a - 2<|> ) ) (6 8) 

(V m +K tj> 2 d )(sin 2a -sm 2$ - 2a + 2<J> ) (6 9) 
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Linearizing (6 8) and (6.9) and , 

v cma arranging in matrix form, we get 


AV a 


' 

CM 

L^_ 

% 

« [ Ql 

^ 2 q 

w ^ 


— w 

•where 



il Aot + i 2 


(6 10 ) 


Q 


= 4 


3IC 
2 4i 


Si =l/| 


3 

St 


4 3 


^2 ® ^ *2^jt 


and 


-VI 


3 


cos 2a - cos 20 sin 2a -sm 20 + 2a - 20 

sm 2a -sm 20 - a + 20 - ( 00s 2 a - cos 2<|) ) 

-( V m +K ^ 2d a + K 0 2q (cos 2a + 1 ) 

K < l> 2q sin 2 a + (? m+ K 0 2d )(cos 2 a -1 ) 

(V m +K ^2d )sin 2<t> ~ K ^2q (cos 2 *+ 1) 

-K 0 2q sm 20 + (V m +K^ 2d )(-cos 20 + 1) 


cos 2 a - cos 2 <J> 

sm 2 a- sm 20- 2 a + 20 + 4n/3 


The term A0 in, equation (6 10) can he obtained from the nonlinear 
equation defining 0 The expression for a 0 as derived in 
Appendix J is 

A0 gd I 

+ g 5 Aa + g 4 A V m (611) 


A 0 =S j £ r ' 


T 


A0 


2q 


jg» g^ and are defined in Appendix J 


.1 
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Eliminating A<f> from (6 10), we get 


+ h,j Aa + hg A V. 


— 


— 

4 v a 

= CHJ 

4* 2 a 

4V q 


_ 4 *2q 


m 


(6 12 ) 


where 




{HJ ssfQJ + t 0,2 &. 1 * H-) ” Si + S2 S3 > hg * S3 *** S2 ^4 

Substituting, equation (6 12) m (6 7), 


! 

0 

H 


r ~~ 

Ato 

2d 


= s 2 


A i„ 

iq 


*2q 


t S-j Aa + Jig ^ 


m 


(6 13) 


where 


fS 2 1 = CS 1 ] + (V 1 H] , ii = 1 ^ li n and 2 2 ^ V 1^2 
Finally, equation (6 5) becomes of the form 

x^Ax + bu + cv+dw (6 H) 

where 

X T = { Ai); 2d W 2q AS| , u = Aa, v = AT m , w = AV^ , 

A = fjL,! + {A 2 } , [A 2 ]= | Cl , b ” * B ll JLl * and 

d = I B 1 1i 2 

let the output vector be defined by 


y = A S = I 0 0 1] x 


(6 15 ) 


f 

A 

A 


\ * 
i r 


i 
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During steady state, the solution of x is given hy 

X(s) = (si- A) + cv(s) + d v/ ( s ) ) (6 16) 

and output Y(s) is given hy 
Y(s) * tO 0 1 1 X(s) 

* u (s) + & 2 (s) v (s ) + Gj(s) w(b) (6 17) 

6 3 AS EXAMPLE 


The A matrix and its eigenvalues and vectors h, c and d 
of the test machine are numerically computed for an operating 
point defined hy a * 70°, S * 0.03 aod V m = 325 V. Load torque 
is considered to he proportional to speed, m this example The 
effect of such a load is to increase the damping coefficient (D) 
from D to (D + T-VoO The computed results are 

f — hi r — i ^ 



-40 868 

-46 200 

103, 

025 

A » 

37 206 

5 845 

-229 

317 


w# 

4 290 

- 2,222 

- 2 

903 



^0 


-0 138* 


o « 


0 

f d * 

0 088 




0 796 


-0 02?_ 



20,92 
-26.11 
6. 30_ 


The eigenvalues are 


X 1>2 = -4.889 ± J 54 87, * 3 = -28 149 

The open-loop transfer functions G-j.Gv, ^3 are 
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, s 6 515 (s^+50s+85l ) 

T1 (s+28 15) (s^+9 78S+1239 8) 

/ v = 0 796 (s 2 +35s+14S0) 

J2 S ~ (s+28 15 )(s 2 +9 78^1239 8) 

& ( s \ _ -0 027 (s 2 +20 48S+828 9) 

5 S ~ (s+28 15) (8^+9 78S+1239 8) 

The block diagram of the motor is given m Pig 6 1 
6 4 STABILITY ANALYSIS 

The stability of the drive under dynamic operation is 
determined from the location of the roots of the characteristic 
equation or eigenvalues of the system The eigenvalues of A 
matrix defined m (6 14 ) are obtained numerically for different 
operating points The loci of complex eigenvalues of the test 
machine (Tabic 2 1 ) arc given in Pig 6 2 for constant value- 
of a It is interesting to observe from Pig 6 2 that there 
is a possibility of instability at speeds very near to rated 
synchronous sneed This uhenorena is rather unexuected and is 
not observed on machines supulxcd by constant voltage source at 
rated frequency [5 7] In order to stuc’y this ohenomena m more 
detail and verify its existence nine more machines are studied 
v/hose parameters arc taken from reference [40] These machines 

can he divided into two groups (1) with (R/V > 1 

with (E g /R r ) < 1 where (R.M r ) is the ratio of stator to rotor 

resistance The parameters of these maenmes are given m 
Table 6 1. Pig 6 3 gives the loci of eigenvalues for the 


Fig 6 1 Block diagram representation of an 
induction motor with voltage control 
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Flg*6.2 Locus of eigenvalues of the test machine 
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Fig 6 3 Locus of eigenvalues of 
machines with o (' >1 




TA^LE 6 1 


MOTOR PARAMETERS 


SI 

Fo 

M* «•# 

1 

1 f-d 

V 

LL 

Poles 

Ereq 

(p U 

)(pO 

K 

m 

(p u ) 


J 

a 

1 

15 

230 

4 

60 

0199 

1270 

2 

073 

0 

150 

1 

489 

2 

30 

325 

6 

50 

0250 

0873 

3 

877 

2 

736 

1 

067 

5 

110 

364 

4 

50 

0170 

1215 

4 

140 

5 

000 

1 

237 

4 

250 

2300 

4 

60 

0141 

0864 

3 

026 

3 

459 

1 

698 

5 

3 

220 

4 

60 

0377 

0349 

1 

208 

0 

044 

0 

533 

6 

3 

208 

6 

60 

0324 

0566 

1 

026 

0 

312 

0 

645 

7 

3.75 

400 

4 

50 

0799 

0526 

1 

127 

0 

1023 

0 

415 

8 

25 

460 

4 

60 

0472 

0498 

1 

951 

0 

277 

0, 

.465 

9 

50 

460 

4 

60 

0402 

0532 

2 306 

0 

331 

0 

382 


2 

Per unit quantities calculated from Z s = Vj,j / (bp x 746) 
All motors are assumed to be wye connected. 
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machine -with (Rg/Ry) * 1 Instability is observed in all these 
cases It is to be observed that the test machine also has 
<VV > ^ ^ 0G1 °f complex eigenvalues for machines with 

(Rg/'ftj.) * "1 ar ® plotted m Fig 6 4 It can be observed that 

the characteristics are Quite different from that given m Fig 
6 3 There is no instability for voltage controlled motors having 
(R g /R r ) less than unity In both cases, real eigenvalue shifts 
to the right half plane of the s-plane when the operating point 
shifts to the negative slope region of torque-slip characteri- 
stics From these results, it may be conjuncture! that voltage 
controlled induction motors with (R /R ) greater than unity may 
experience instability very close to the synchronous speed 
However, this phenomena may noo be observed m practice since it 
occurs for a small region near synchronous speed 

loci of the complex eigenvalues of the test machine with 
sinusoidal excitation obtained from a linearized third order 
model (Appendix K) are shown m Fig 6 2 for different input 
voltages. For a voltage controlled induction motor with sinusoi- 
dal excitation, there is no unstable region as can be seen from 
Fig 6 2, although the value (R g /R r ) is greater than unity. The 
complex pair of eigenvalues always lie in the left half of s-plane 
The dissimilarity of characteristics between the sinusoidal exci- 
ted operation and thyristor fed, voltage controlled operation 
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Fig 6 4 Locus of eigenvalues for machines 
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can be ascribed to the differences m the two cases given 
below 
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1) The effective terminal voltage in. the case of voltage con- 
trolled induction motor is a function of the rotor flux 
linkages while for a sinusoidally excited motor it is 
independent of the rotor variables This is noticed from 
the fact that Av d and 6v q are functions of A'J'ga 311(1 A ^2q 
(see equation 6 10) 

2) The variation of the rotor flux linkages with slip follows 
a different pattern in the case of thyristor fed motor (see 
fig 5.5 of Chapter 5) 

3) The effective terminal voltage of a voltage controlled motor 
changes rapidly with small changes in the slip when the 
motor is operating near synchronous speed (see fig 5 6). 

Since the region of instability is very small for the test 
machine, this phenomena could not be observed experimentally. 

6.5 DESIGN Of THE CLOSED-LOOP SYSTEM FOR SPEED CONTROL 

For precise speed control applications, a closed-loop 
system with speed feedback is essential In this section, a 
voltage controlled induction motor m a closed-loop configuration 
with speed feedback is discussed The transfer functions and 
block diagram of the open-loop system developed in Sections 6 2 
and 6.3 are used for the design of the closed-loop system 
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6 5.1 Description of the System 

A squirrel cage induction motor (test machine) is used as 
a drive with a mechanically coupled d.c generator for loading, 

A tacho-generator coupled to the induction motor is used for 
obtaining a voltage signal proportional to speed The output 
of the tacho-generator is stepped down by means of a potential 
divider arrangement The output of the tacho-generator has high 
ripple content due to the commutator action and it is filtered 
by means of a first order R-C filter This filtered output 
which is proportional to the shaft speed is compared with a 
reference voltage corresponding to the set speed of the drive 
The error signal is fed to a proportional plus integral (PI) 
Controller to adjust the triggering angle a such that the error 
signal at the input of the PI controller tends to become zero 
A precision half wave limiter 1511 is used at the output side of 
the PI controller so that a always lies in the normal working 
range (0 to 180°) The block diagram of the closed-loop system 
is given in Pig 6 5 

PI Controller for the Speed-loop 


The functional diagram of a PI controller is given in 
Pig 6*6 and the transfer function is given by 


knO+sO 

G 0 ( S ) = — g — 


(6 18 ) 


where 


T c = °o E c2 ’ 


'c cl 
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A half wave preoision limiter -with unity gain is used at the 
output side of the PI controller to limit the voltage within 
the allowable limits, 

Tacho- generator 

The tacho-generator output is stepped down to 3V using a 
potential divider arrangement corresponding to the maximum 
speed of the motor (S * 0) The output voltage is assumed to 
vary linearly and the transfer function is given by 


*«<■> - -43 - - 3 (6 195 

Filter 

A first order R-C filter (Pig 6 7) is used to reduce the 
ripple content of the tacho-generator voltage The transfer 
function is given by 


s f (s) " i+sT 


(6 20 ) 


where 


T f - R f 0 f 


Triggering Circuitry 

The triggering angle a vanes between * and 0 corres- 
ponding to an input voltage variation of 0 and 57 and hence 
the transfer function is given by 

a ( Q ) * MkL-T " */3 * -1 

Ae (sj ' 


4 


(6 21 ) 


6 5.2 Controller Design 


let the operating point be defined by a * 70° and slip 
S = 0 05 The numerical values of the motor transfer functions 
are given m Section 6 3 for the same a and S The design of 
the closed-loop controller involves the determination of gam 
and time constant of the PI controller and filter used m the 
feedback path 

The filter time constant T^ is kept large enough to reduce 
the ripple amplitude and small enough to keep the cut-off fre- 
quency 1/T^ greater than the gam crossover frequency a> c . 

Accordingly T f is chosen as 8.33 msec 

Prom stability consider at 10 ns, the Pl-cont roller time constant 
T q is chosen such that the gam crossover frequency occurs at \ 

-1 slope and 1/T c < « c < 1/T f C 58,591 Accordingly T c 
is chosen as 250 msec 

The open-loop transfer function of the complete system 

0 45 k c (1 + f)(1 + -2§-3 + gfy) 


becomes 




t 
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Let 


GrH(s) = L'-1 (s)) B ( s ) 


where 


k’ * 0 45 k. 


k 


12 5 


jr 

j 


A(s) ~ 2} » B ( s ) - 0 + s / 2s 5 + b 2 /851) 

0(b) - (1 + s/130 + s 2 / 1185 6) 

Bode plots of A(s), B(s), C(s) and GH(s)/k' are shown m Pig 6 8, 

The gain crossover frequency for a phase margin of 50° is 60 rad/ sec 
and the value of k' corresponding to this value is 15 dB There- 
fore, 


6 6 DETEi 1INAII0N OE DYNAMIC RESPONSE OP THE CLOSED-LOOP SYSTEM 
661 Experimental Results 

On the Basis of the design described in the previous 
section, an electronic circuitry was developed for the controller 
(see Pig 6 9) and experiments were conducted to determine the 
dynamic response under closed-loop operation The Pi-controller 
gain is chosen as 6 although the designed value is 12,5 Mai- 
operation of the triggering circuit is observed for gam greater 
than 6 The inability of the circuit to respond to fast input 
variations may be the cause of mal-operation Photographs of the 
dynamic responses due to load perturbations and change in 
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reference voltage are given in Pigs. 6,10 to 6,12, Test results 
for both positive and negative step inputs m load torque shown 
in Pigs. 6, 10 and 6.11 indicate that the controller is effective 
during both acceleration and retardation of the motor. It can 
also be observed that the steady state error in speed due to 
perturbations in load torque is practically nil. When triggering 
angle ot becomes less than or equal to <J> , the controller be- 
comes ineffective because any further advancement of ot does not 
have any effect on the motor voltages. In Pig. 6.12, the res- 
ponse due to changes m the reference voltage is given. The res- 
ponse is well damped and the speed settles down to a new value 
with little overshoot. 

In Pig. 6. 13, the photograph of the experimental set up is 
shown. 

6.6,2 Analytical Results 

Transient response of the closed-loop system for perturba- 
tions in load torque and reference voltage is obtained by solving 
the system differential equations. The closed-loop transfer 
function is transformed into a set of first order differential 
equations. These equations are numerically solved using fourth 
order Runge-Kutta method for a unit step input of the disturbance. 

The open-loop transfer functions given in Section 6.3 for an 
operating point of S = 0.05 and a * 70° are used for the analysis. 
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Fig. 6.10 Response due to load perturbation (positive step 
input) (25 rpm/division, 2 sec/div ) 




Mg. 6.11 Response due to load perturbation (negative step 
S input) (25 rpm/division, 2 sec/div) 



Por the closed -loop system shown in Pig. 6.5, the closed- 
loop transfer function is given by 


x 

u 


(s) 


V 5 + + + + h A s + b £ 
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~T 


s + ai s + a 2 s + a 5 s + a 4 s + ^ 


The equation (6.23) can he transformed into state 
space form [60 ,61 ] and is given by 
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where c’s are defined as 


C o * b o ; c i = \ ~ 


i-1 

l 

m=o 


a. c 
l-m m 


(6.23) 


u 


( 6.24) 


and x is given by 

x = x-i + c o u (6.25) 

Response (aS) due to Perturbation in Reference Input 

x and u in equation (6.23) become AS and Ae ref respecti- 
vely during this case. The coefficients of the variables for 
the operating point of a = 70 and slip = 0.05 are 
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a1 = 157.9 

b Q . 0 

a 2 = 9476.9 

b 1 = 9.47 

= 552685.0 

b 2 = 1458.4 

a 4 = 7499927.0 

b 3 = 47825.5 

35 = 11607980 

b 4 = 115568 


b 5 = 5868505 


The solution of equation ( 6 . 24 ) is obtained for u = 1 ■with all 
initial values set to zero. The computed result ( AS) is shown 
in Fig. 6.14. 

Response (As) due to Perturbation in A T 

■ — > — — m 

The coefficients m equation (6.23) can be obtained from 
the closed-loop transfer function (see Fig. 6.5) relating S to 
T m and are given by 


a 1 - 

157.9 

b . 0.0 

0 

a 2 = 

9476.9 

b 1 = 0.796 

8.-2J = 

552685 

b 2 = 125.4 

a 4 = 

7499927 

b 3 = 4521.5 

= 

11607980 

b 4 = 141569 


For a unit step input inAT^, the solution of equation (6.24) is 
obtained numerically. The variation of As with respect to time, 
as computed, is given in Fig, 6.15. 



Time (sec 


Fig. 6.14 Response due to refrences speed perturbation 



Time (sec 


* i i 

Fig. 6 15 Response due to load perturbation 




6.6.3 Discussion 


1 68 


The solution of equation (6.24) is obtained numerically 
by fourth order Runge-Kutta method. A step size of 0.0 2 sec. 
is used for the integration. The computed result (Pig. 6.14) 
for a disturbance in A e re f shows that AS reaches a steady 
value of 0.33 in about 2,2 seconds. Comparing with the experi- 
mental result given m Pig. 6.12, it can be observed that al- 
though the responses are not alike the time needed to reach 
steady state is nearly same. 


The response (AS) due to unit step input in A T m , as given 
in Pig. 6.15, indicates that there is an initial change in slip 
of 0.02. The system gradually settles down to the initial 
value of slip in 2.5 seconds. The results of experimental stu- 
dies (Pigs. 6.10 and 6.11) show good correlation with the ana- 
lytical results. The system could not be operated at the 
designed value of the controller gam and the operating value 
had to be reduced. This may be because the dynamics of the 
triggering circuit is not taken into account when the controller 

gain is designed. 


6.7 CONCLUSIONS 

A third order, small signal model of a voltage controlled 
induction motor drive is developed. It is observed from the 
eigenvalues of the system that instability may occur near about 
maximum speed for motors having ( R s / R r ) greater than unity. The 
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CONCLUSION 
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7.1 GENERAL 

The major aim of this thesis has been to provide a 
comprehensive analysis of the steady state and transient 
performance of induction motors with voltage control us in s 
thyristors. The application of voltage controlled induction 
motor drives requires a proper design and selection of the 
motor for reliable operation and hence, the knowledge of the 
transient and steady state performance of the drive is of 
great importance. The earlier approaches to the steady 
state analysis are based on state space techniques and har- 
monic equivalent circuits and the limitation of these appro- 
aches is either the lack of generality or the computational 
complexity. Very few articles dealing with the analysis 
of transient performance of voltage controlled induction 
motor drives during open-loop and closed loop operations have 
appeared in literature. This is mainly because of the lack 
of suitable models for the analysis. This thesis has presen- 
ted two novel hybrid models of an induction motor. One of 
which is well suited for the transient and steady state 
analyses of voltage-controlled induction motor drive under 
wide range of operating conditions. Specific contributions 
made in the thesis are reviewed in subsequent sections. 
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7.2 STEADY STATE ANALYSIS OP THYRISTOR 7CLTASZ SORT ROLLED 
DRIVES 

The generality of the B7A in the steady state analysis of 
voltage controlled induction motor drives has been illustrated 
in Chapters 2 and 3 where the steady state analysis and simu- 
lation of a motor with 

a) a pair of "back-to-back connected thyristors in each 
line (Chapter 2), 

b) delta connected thyristor bank at the open neutral 
of the stator winding (Chapter 3), and 

c) back-to-back connected thyristor-riiode pairs in each 
line (Chapter 3) 

have been presented. In the BVA approach, the motor equations 
are invariant while the different inodes of operation, are 
accounted by altering the forcing voltages. The open 
circuit condition is simulated by applying an equal but 
opposite voltage to the open circuited winding so as to 
force a current zero condition in the winding. 

A comparative study of the performance of the drives 
with different thyristor configurations has also been 
presented. 

7.3 CIRCUIT MODELS OP INDUCTION MOTORS 

Por easier formulation of thyristor controlled drives, 
it is necessary to have a three phase circuit model of the 
machine, which retains the identity of the motor teminals 


4 
% , 
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and yet is computationally simple. Two novel hybrid models 
of an induction motor have been described, one particularly 
suitable for stator controlled machines vhile the other is 
well suited for rotor controlled machines. Simple, three 
phase equivalent circuits have been presented which can 
accomodate any external terminal constraints. The need for 
the inversion of a time dependent matrix is eliminated in 
this model. Using one of the hybrid circuit models, the 
formulation of state equations for different terminal condi- 
tions has been described with the application of network 
theory . The simulation results indicate that the rotor flux 
linkages referred to a synchronously revolving reference 
frame are essentially constants during steady state even when 
the stator current contains large harmonic components. This 
implies that the constant component of the torque , as per 
the model presented in Chapter 4, is produced solely by the 
fundamental component of current during steady state. This 
fact is of great significance m developing models suitable 
for the analysis of closed loop systems. 

7.4 PREDICTION OP STEADY ST APE CHiRACT ERISTICS 

Determination of the steady state characteristics of a 
voltage controlled induction motor drive generally involves 
numerical integration or evaluation of state transition 
matrix over a certain duration. In Chapter 5, a simplified 
analysis has been described by which the steady state chara- 
cteristics such as average torque, phase angle 4> , fundamental 



173 


components of stator current and voltage can be obtained 
accurately by solving tnree nonlinear algebraic equations 
by Newton-Raphson method. Computationally, this method is 
much faster compared to other methods* 

Another advantage is that a direct relationship between 
the motor performance and control variables is available. 
This is necessary for the analysis and design of closed-loop 
systems for speed control applications. 

7.5 STABILITY ANALYSIS 





No work has been reported in the previous literature on 
the stability aspect of voltage controlled induction motors 
using thyristors. Prom eigenvalue analysis. It is observed 
that the induction motor drive experiences instability dur in g 
a narrow range of operation close to synchronous speed. Prom 
a case study of 9 machines, it is observed that machines with 

(Rg/Rj.) greater than unity show t tendency for instability at near 
maximum speed when subjected to voltage control using thyri- 
stors. 

7.6 DYNAMIC ANALYSIS OP THE DRIVE WITH PEEDBACK 

Using a linear model of the drive and its associated 
transfer functions, a closed-loop system with speed feedback 
has been designed and developed. Dynamic response due to 
perturbations in load torque and reference speed has been 
obtained both analytically and experimentally. 


Ph • 
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7.7 SUGGESTION EOR FUTURE WORE 

®he following topics are suggested for fut\^ re ^ 

1) One of the hybrid models developed in Uhapt ev , 4 has 

been used for the digital simulation of v oltage oontrollea 
induction motor drives. The same model >, e use ^ f 0T 
the digital simulation of other types of st a t or contro- 
lled machines such as inverter-fed drives. 

Although a model suitable for rot or con t ro iiea 
machines is developed in Chapter 4# it has riot been uged 
since it does not come under the purview 0;f thesis 

This model can he used for the digital s imux a ti 0 n of rotor 
controlled machines. 

2) A detailed study of the stability aspect 0 f voltage 
controlled drives may be undertaken to inve s txg a t e the 
causes of instability observed in certain category of 
machines. 

3) The study of the closed-loop system presented in 
Chapter 6 is of preliminary nature. A comp rebeasive 
study of the closed-loop system Incorporating 0 ther 
types of feedbacks such as current feedback ari .^ a 
detailed analysis using modem control theory can he 
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JG?j? HiJD £ A 


boundary RELATIONSHIP BETWEEN miAL LD INITIAL 
VECTORS AT A DISTANCE OP if / 3 


Since the voltages v ag , and v cg are nal- 

/-mmetric 


v as^ t + ^ = “ v as 


(«i»t + * ) = ^ 
DS 


Ids 


v cs < ^ t + ir) = ~ Vcs 


U.1) 
(A. 2) 
(A. 3) 


addition, these three-phase voltages are displaced by 
r j 3 radians. Hence, 

(A*4) 
(A. 5) 
(A.6) 


v bs (<rt + 2V3) = t as (»t) 

v (lot + 2v/3) - v te 0»t) 
cs 

v (<rt + 21I /3) - v C»t) 

as 


creasing the arguments of equations (A.4 - *.6) W ^ 
a substituting the results in (A.1 - A.3) y ields 

v as ( B t + «/3) - - v bs (,,rt) (A ' 7) 

v (lot + »/3) * - v os^“ t ^ 
bS (A. 9) 


, (ut + »/3) = -’as 0410 

CS 


, „ ^ a and phase variables are 

relationship between 4-9 an4 ? na 



1S2 

'as *= 

(A.10) 

- - 1//6 V ds " 1//2 

U.ii) 

Y os = -V/6 v ds ♦ V/2v qs 

(A.12) 


Substituting equations (A.10 - A.12) in (A.7 - A. 9), we get 


V ds 
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✓3/2 

T" 1 

P* 

Vqs > 
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(A.15) 
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Similar relationship holds good ior the variables Y s1 , V sZ , 


v r1 3114 v r2* 


1. e. , 


(Me ^ 
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7 r1 


z*. 

ft.*. Jkm. 
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•SEPEBD II B 


BOUNDARY SEMIOSSHIP BHfM® TS3 FINAL AND 
INITIAL VECTORS AI A EIS? ARCS OF 2ir/3 

Ihree phase symmetry demands that 
v bs («t + 2«/3) = v (mt) 

SIS 

v 03 4.t + 2»/3) = T bs («t) 

v as (“t + 21T/3) . (ut) 

cs 

The phase voltages are related to the d-q voltages by 
"as “ ^ 

v bs - -V/6 v dg . 1//2 v qs 

V OS * V 4g + 1//2 y qs 

Substituting equations (B.4 - B.6) in (B .1 - B.3), we get 
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( 3 . 2 ) 
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Similar relationship holds good for the vector f 7 „ V v V 1 T 
rr. t i— ^ si s2 rl r2 J 
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APPENDIX 0 


ELECTROMAGNETIC TORQUE IN A MOTOR WITH 3-PHASE 
STATOR AND 2-PHASE ROTOR 


The general expression for the electromagnetic torque i£ 


T = 1/2 i T i 

e 36 r 

Referring to equations (4.1 - 4.4) 


(C.1) 


m m rp rp 

1 = *V VI 5 x s = 1 i 1 a i W » 


I o n ; 


T . . . 3 L 

x r =ll 2a X 2 6 J ; 39, 


T-sin 0. 


cos 0. 


/ oj 

1 

* 


L 12 -sin(0 r -2Tr/3) cos (0 r -2w/3) 

-sin(9 +2^/3) cos(0 +2ir/3) 

i_ r j. „ 


Therefore. 


T - il i 1 i T 1 I 

x e “ ^ 1 1 s r 1 | 


f 0 5' 


T e - 1 3 I h i r 


Equations (4.18) and (4.4) are rewritten as 
i s " °1 i 1dq 

V 


(C.2) 


( 0 . 3 ) 


( 0 . 4 ) 



1S5 


■where 


rrrb 


it* 

cos 9 q 

cos(9 0 -2lf/3) 

cos(6 q +2V3) 


sin 9^ 
o 

sin(0 o -2n/3) 

sin(9 Q +2V3) 


.? t ,• i 1 

* 1 1dq = 1 x 1d IQ 1 


T T 

vr * i» 12 


cos 0 r cos(0 r -2V3) cos(6 r +2V3) 

sin 0 r sin(6 r -2^/3) sin(0 r +2*/3) 


*2$1 


The relationship between o - 6 and 1-9 variables is 


(0.5) 


♦ r - G 

•where 


2*2dq 


cos(9 o - 0 r ) sin ^o” 6 r^ 
-sin(9 0 -e r ) cos(0 o - e r ) 


to 

* w 2dq 


2d ’ J, 2q 1 


, c S) in (0.2) and after simpli- 
Subst itutmg equations (0.^ - .5) 

fication, we get 


* e -/?7* 


1 2 / • $ 
17 U 1d ^ 


- ij lq *2d) 


( 0 . 6 ) 
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APPENDIX D 

RELATIONSHIP BETWEEN THE PARAMETERS OP THE HEBRIE 
MOTEL MD COFTENTIONAL PjSOTE EQLT/.'iEII CIRCUIT 


Let the two-axis rotor currents in the hybrid model be 


and 


, .r 


cos (e r - e # ) 

(D.1) 

r 5in ( 0 r - 9 o ) 

(33.2) 


The stator flux linkage is given by 
♦ | = tt 1 -l 11 )l 1 a +I'12 <=° s * r X r Qos( V V* 

1 12 sin 6 r I r sin(e r - 9 0 ) 

- - ln^la ♦ x 12 r r 008 0 o 

= (I-) - L 11^ii a + l 12 *-0 

= (L 1 - i-liU-ta + 1*12 

where i r and i* are the rotor currents referred to the stator. 
Referring to the conventional passive equivalent circuit 
(Rig. D.1(a))» the stator flux linkage is 

s • ($• 4) 

n> s a- L i 4 + m l 
a s » a a 


(3). 3) 


■where 
m = 


, L g = X s /« * ^ + X dL )/<0 
Comparing equations CD. 3) eJld CD. 4) 





1S8 


L 1“ L 11 3 L s 

It is possible to prove 128,62] that 



(2.5) 

( 2 . 6 ) 


(2.7) 


Therefore, 

2 w (2.3 

l i - % \ * x a. )/u 

Prom the definition of L r (equation 4.4), it is evident that 

L 2 - Ij/* « ( x m + x rl^“ 

L s in equation (4.11) is defined as 
' 2 

3*12 

L!f - h - L 11 - Tip 

> « s - sX )/u 


(2.10) 


189 


APPENDIX E 

TORQUE EXPRESSION EOR A MOTOR WITH 2-PHASS STATOR 

AND 3-PHASE ROTOR 

The general expression for electromagnetic torque is 


T = 1/2 i T i 

e 3 8 r 


Referring to equations (4.22 - 4.26) 


(B.1) 


i = lig 1 » ig * [ 0 *i g 1 » 1- r = f *2a ^2b i 2c * 


3L 

30, 


0 1ST* 

u ,T 0 


, JP-* Li 


12 


-sine r -sin(6 r +2tr/3) -sin (6 r ~2it/3) 
-cos 8„ -cos(e +2»/3) -cos(0 -2*/3) 


Therefore, torque expression becomes 


T 

T = i x N i_ 
e s r 


(E.2) 


The relationship between the phase currents (rotor) and d-q 
currents is 


i r “ C 3 i 2dq 


(E.3) 



C 3 =/273 


cos(e 0 -e r ) 


sin( 0 Q -e r ) 


i* 

* 1 2dq 


* i 2d *2q* 


where ^ 

cos(0 o -0 r -2»/3) sin(8 a -ft r -2ir/3) 

cos(0 o -0 r +2w/3) sm(8 o -0 r -i 

m 

Equation (4.27) defining the stator current i g is rewritten as 
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where 


M' 


JE 


L i2 


■s ’ rf<* B - mT l x > 


cos 9 r cos(e r + 2 ir/ 3 ) cos( 0 r ~ 2 it/ 3 ) 

•smfl r -sin( 9 r + 27 r/ 3 ) -sin (® r - 2 ^/ 3 ) 


(E.4) 


+s 1 » l * 1 « ♦in’ ’ 


The relationship between a - (3 and d-q stator flux linkages are 


*3 - °4 +iaq 


(B. 5 ) 


where 


cos 6 sin 0 ( 

•sin 9 _ cos 0 ( 


^Idq =t ^ Id 


_ D 

Substituting equations (E .3 - E. 5 ) in (E. 2 ) we get 

t = m 1,2 
0 


*14 i 2q " ^Iq X 2d 


) 


(E.6) 
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APPENDIX E 


DERIVATION OE STATOR CURRE1\[T EXPRESSION 


Eor the purpose of analysis, the half cycle period of 
the current waveform is divided into six sections as shown 
in Eig. 5.1. Referring to equations (5.1 - 5.4), the current 
expressions for each section are derived as follows : 

Section 1 a < ® 0 i (60°+ 4>) 

The differential equation defining the stator current is 


Rl 1- + 1? p(i 1a + I ia )- \ sin9 o 


(E.1 ) 


where 


6 = wt 

o 


Rearranging the terms and replacing pl^ a (as defined in 
equation ( 5 . 2 )), 

Li* T>i+ + L i-i * E sin(0 + e ) 

1 y la 1 la 0 


(E.2) 


where 

E = (Eq + , E q = V m + i “d 

w I 


K 


-2 ^ p2^^ 

J 12 


+ s * 2a . E a = ' K *2q' 
t .ii 

0 s* tan jj 


Solution of (R.2) with the 


initial value of i 1a - 0 at 


0 = a is 

o 
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i 1a ** |-frrj- { sin(0 o + 5 ) + 0 1 exp(-B(9 Q -a) )} (M) 

where 

R-* -1 1 

Zl| = H| + 3 X|, 6 = 8 - p, 0 = jjr , P = va n (— ) 

X^’ =wLl| and 0^ = -sin (a +5) 

Section 2 (60°+ 4>) <, ® 0 1 (60° + a ) 

In this section the current in phase <5 is zero. The 
initial value of i-i during this section is the final value 
of i, in the previous section. is i.e., 

i 3- 

I 2 = ^|np{sin(ir/3 +<J> + 0 + G 1 exp(-0(tr/3 + *-<*)) > (1.4) 

Differential equation defining i 1a is 

E, l 1a + L!j P i 1a - /V2 B state q + »/6 + 9) (?-5) 

On solving this equation with I 2 as initial current, we get 
! 1a = -j§,|{/3/2 9i“( 0 o + l + {) + °2 eipC-BCV' 73- * ))>(I '- 6) 
where 

0 2 . 0 1 «cp(-K»/5+* -■»)) + +*+«)- •'V 2 cosC^+6) 

Section 3 (60° + «) * e 0 - (120 ° + * ) 

In this period all the 3-phases carry current as in 
Section 1. The differential equation defining i 1a is same as 
(p.2) with an initial current given by 
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I ^ = i ' Tg rr j{ / 3 / 2 cos(a+ 6 ) + 0 2 exp (- $ ( a - <p )) ) *7/ 

The solution of i 1a is given by 
a t= y^T Tj { sin (9 exp (- 6 (9 Q - ^ /3 - ®) ( x * 

■where 


C 3 * C 2 exp (- 0 (a- 4> ) ) + /3/2 cos (a + 6) - sin(ir/3 + a + 5 ) 
Section 4 (120° + <!>)< 9 0 < (120°+ a) 

Phases a and c are conducting while h is open during 
this section. The initial values of i 1a is obtained from 

2 IT 

(P.8) by substituting e Q = y + * • 1 » e *» 

I 4 » y|y| { sin(2ir/3 +4+ 5) -r C3 exp(- 9 (»/3 + * >> (3? * 9) 

The differential equation defining i 1a is 


L i' P ba + H 1 


i ./3/2Esinfo - ir /6 + e ) 
la 0 


(p.10) 


Solution of (I'*10) is 


jL _ -jL.{/ 5 / 2 sin (9 - ir/6+ «) + O4 exp(-9 ( 9 0 - 2ir /3 - 
la }2^j (P.11) 


where 


C 4 . C, expC- «(-/3 ♦ ^ -“<^ 3 + * + S) - /3/2 C0S( * +6) 

Section 5 0 2 °° + ° 5 - ®o * l1 ®°° + * 5 

Ita ee phase conduction tafces place as in Section 1 and 5. 
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The initial value of current during this section as obtained 
from (F.ll) is 

I 5 = yfirj f ^3/2 cos (a + « ) + exp(- p ( a - 4 > )) } (3 s . 12) 

Since the differential equation defining i^ is same as (3 s . 2), 

*21 

solution of i .4 is 
» a 

i 1a = { sm (8 0 + <5 ) + exp(- 6(8 Q -2V3 - a))} (?.13) 

•where 

exp(- 0(a-$)) + /3/2 cos(<M-d ) - sin(2ir/3 + a + O 
Section 6 (180° + 4> ) < 0 Q < (ISO® + a ) 

During this period i^ a is zero since a-phase is open. 
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APPENDIX G 

CALCULATION OE FOURIER COMPONENTS OF THE VOLTAGE WAVEFORM 


The general form of Fourier series is 


- f( 9 0 ) = 

a + 
0 

00 

l (a n cos n e Q + b n sin n 9 0 ) 



n=1 

where 



1 

a o = 'Si 

2 ir 

J 

f(e o ) a e 0 


2 V 


-1- 

I! 

1 

0 

f (0 0 ) cos(n 9 0 )d 6 0 

-r 

2 ir 

/ 

0 

f(8 0 ) sin(n 9 O )^0 O 


The voltage wave shape of Fig. 5.1 has zero time average 
value so that a Q Is sen, In ( 6 . 1 ). The fundamental component 

a-] is given by 
(ir + a) 

1-f I 7 as “ S 0 o d e o 

a 

Let 

a 11 = fundamental components of voltages In Sections 1,3 and 5 

a, 2 = fundamental component of voltage in Section 2 

a,, = fundamental component of voltage in Section 4 

a 14 = fundamental component of voltage in Section 6 

so that 

a 1 - a i1 + a 12 + a 15 + ^4 
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a 11 = 


(ir/3 + $) (2 tr/5 + <fr) (ir + 4 ) 

$ V m aln «o 008 a 9 o 

a 
a 


(w/3+a) (2*/3+a) 


= 0 


{C 2) 


(n / 3 + ia$ 


a 1 2 * Jf ^ (0 o +ff/6 ) + ®i smfe 0 +2ir/3) Eg cos (0 o +2tr/3) } 


(ir /3 + <}> ) 


cos e o <n o 


= ^ { l cos(2a+2ir/3) cos(2<t+2ir/3) J [ 3 V m +B^ /3 E 2 ] + 

(sin(2a+2ir/3) sm(2<|> + 2if/3)l {/ 3 V m+ /3 B 1 + E 2 ] + 

2 (a «) 1/3 V m + /3 E 1 + E 2 1 > (G 3) 

(2if/3+a) 

= 2 f{ n V Slate it/6)^ sin(0 o 2*/5) B 2 008(6^211/5)) 


a 1 3 53 F2 J ' ' ' m o 

(2ir/3 + <t*) 


cos e 0 d e 0 


a 


. 1— r loos (2 <*h 4 tr/ 3) cos(24»+4‘ ,r /3)l t + ®1 + 3 

8 ir 

l sin (2t* +4*/ 3) - sln(2$+4*/3)1 l ^3 v ra /3 B 1 + E 2 1 

2(a-6) I /3V, /3 3, + % H (S4) 

(n+a) 

1 J { sin e 0 + 008 e o* 009 0 ° d $0 


14 “ it 


(*+<>) 


. y (008(20 ) 008(26)1 E 1+ 3! 2 (8xn(20)-sin(26) + 2(«-*)l^ 
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Henoe 


a 1 = + a 12 + + a-, 4 

» Tpjj, t (cos 2» cos 24 ) (V m +E^ ) + E 2 (sin 2 a sin 2 4> + 2 a - 2$) > 

(G 6) 

The fundamental Eourier coefficient is defined as 
(it +«) 

‘l I v as sln 9 o 4 9 o 

a 

“ b 11 + *12 + b 13 + b 14 
(t/3+«) (2ir/3+*) (tt+4) 

*H - | / V m sin e o 3111 e o d 6 o * 1T {TT 3(a (G 7) 

a (ir/3+a) (2ir/3+«) 


b ^ 2 *° *14 can * e similarly derived as 

*12 “ ■§* ^ f oos(2ci+2’( t /3) cos(24+2tt/3)j t ^ ^ E 1 EgJ + 

(sin(2a+2ir/3 gin(24 + 2ir/3)] t 3 V E r /3 E 2 ) + 

2(a ♦) 13 V m E 1 + /3 E 2 i> < G ! 

b 15 * ( cos ( 2 a +4 V 3 ) - cos (2^+4 it/ 3)1 C /3 Y m +/3 S 1 B 2 J + 

t sitt(2a+4 ir /3) sin (24 +4 1T /3)l l 3 V m +E 1+ i/3 E 2 1 + 
2(o-«>) (3 V m - E 1 / 3 E 2 J > ( & 
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-j 

*14 * 7 ^ E 2 [cos( 2 a)-.cos( 2 <j.)) +E 1 {sin(2a^ sin{2<>) 2(ot-<>)l > 

(G 10) 

Hence 

3 

*1 = Tfi t E 2 [cos(2a) -cos (2<J>) j + 


(W {sin(2a) sin(2«J>) 2 (a 4>)j + ? m } (G 11) 


The coefficient may he expressed as 

3*1 = a 1n + a 2n + a 3n + a 4n + a 5n + a 6n n / 1 

(*/?+♦) 

a 1n = ” J V m sin 6 o 005 n e o d *o 

a 

On simplifying we get 

jn f cos (1 n)(it/ 3 +<fr) cos (1 n)a 
a 1n " n 1 (1 n) r 


cos(1+n)(i|/3+i|i) cos(1+n)a j ( G 12 ) 

(1+n) 

(V3+«) 

a 2n " / {/3 V m sin (e o +1I /6)+E 1 sin (0^2*/ 3) 

(n/3+$) E 2 cos(9 0 +2V3)) cos n 0 Q d 9 0 


^fU -3 V B 1" /3 E 2 ! , a °" (1 t 

oos(1+n) (tt/ 3+<*| ^ ^.cos(1+n) (ir/3+4>lj + 

,/ 3 I ♦ 

flin(1.t.n)(ir/3 *«) sin(Hn)(ir/3+»lii (0 13) 










cos(1+n) (tr+a) oos(1+n) (ir+ 


sin(1+n ) (ir+a) sin(1+n 


3in(1+n) ( tt / 5+ 6) sin(1+n 


{ i # Y m E 1+ /3 E 2 1 t S1P<1 


m(1+n) (tt/ ^+ a) mO+n) (iy3+») - 


t /3 V m+ /3 E 1+ E 2 ] l 00 ^ 1 ““s' 1 * 

cos(1+n) (tt/3+q) cos(1+n)v , w/3 +4> 


Bin(1 n 


sm(1 a)( 1T /.5± 
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APPENDIX h 


RELATIONSHIP BETWEEN ACTIVE AND PASSIVE EQUIVALENT 

CIRCUITS DURING STEADY STATE 


From the circuit model of Eig 5 2 
be written as (see eqn (5 °^) 

the current 1 can 

^ I3X 1 

1 - wn 

7 1 

(H 1) 

The rotor flux linkages as defined in Section ^ of Chapter 4 

are 


^2d = R 2 *2a S ^2 * 2q * R 2 m 1 a 

<H 2) 

V * 2g = s x 2 * 2(J r 2 * 2q + r 2 n lq 

(H 3) 


where 

m » /3/2 I 12 ( see Appendix D) 


Since and are const ant g during steady state the 
terms p and p ^ 2 q can be set equal to zero in (H 2 and 
H 3) The resulting equations are 


R, 



+ S A z 



h " x i 


(H 4) 


R, 



3 X 2 ^2d R 2 m x q 


(H 5) 


Multiplying (H 5) by ;} and adding to (H 4) we get 


* 


*2 m 


(H 6) 


where 


Rg 4* 0 ® ^2 
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Tho relation, hip “between I and ^ as given m equation (5 Vjr) 
is 


r m 

I_ W ♦ 


(H 7) 


Substituting equation (H 6) m (H 7) we get 


X =* 


W7 


(H 8) 


Hence equation (HI) becomes 


i = (H 9) 

w m R 9 

(2 1 + 

Let u now consider the conventional passive equivalent circuit as 
shown m Pig D 1 of App D where D 1(a&b) are equivalent 

The impedance 2^ Z 2 2 and are defined as 


Zfi = + JwL.j 

Z-12 = m 

Z 22 a r* + 3<1) L 2 


where 

L 1 - < X sl + *2 - (x ri + X m )/w m = V“ 


Prom Tig I) 1(b) 


Y sb + Z-|2 -^2 

(H 10) 

0 - E-j + ^22 *2 

(H 11) 
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Eliminating I 2 from equations (H 10 and H 11) 

z 2 

V = M Z 11 Z^ } 

i e 

- T/Zj 

where 


= Z u Z 12^ Z 22 


= Rj + 3 w 



Rearranging the terms after adding and subtracting a term 
(3wm) 2 /owE 2 to the right hand side of the equation we get 

0 2 V s 

Z T « R, + 0 M (L-j rn Z /L 2 ) + (3»o) < j »T£) } 

2 

in Appendix D it is proved that (L 1 ® l^z) ” L 1 

Therefore 


Z T * Z 1 + 3 « ( L 2 Z 2> 


i e 


I, = 


1 Z-j + 3 m rn Rg/ (L 2 Z 2^ 


(H 13) 
(H 14) 


Comparing (H 9) ana (H 14) it ran be seen that i - X, 

Hence the active circuit given in Pig 5 2 is equivalen 
the conventional passive equivalent circuit (Pig D O ««“« 

sbeedy state 
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APPENDIX I 

RELATIONSHIP BETWEEN HARMONIC COMPONENTS OP D Q 
AXIS CURRENTS AND PHASE CURRENTS 


The transformation matrix relating d-q axis and phase 
currents is 


= /f 


cos 8 q cos(e 0 2ir/ 3) cos(0 o +2tt/3) 


sin B q si n(e o 2ir/3) sin(G 0 +2ir/3) 


(I 1) 


For a nth harmonic phase current the d-axis current can he 
expressed as 

i d (t) *s S 2/3 { sm(ne o +c» 11 ) cos0 o +sin(n(9 o 2ir/3)+a n )cos(9 0 2rr/3) 


+ sin(n(8 0 +2V3)+a n )cos(8 0 +2ir/3) > 


(I 2) 


where 

X n = * b n^ a n and h n are Fourier components of 

phnuo current 

Equation (I 2) after simplification becomes 

I d (t) \ I n { [sm((n+1) V a n) + sin( ( n+1 )( e 0 +21T/3) + a n )+ 

+ sm((n + 1)(e 0 2^3)+%)! +{sin((n 1)fl 0 +« n ) + 

sxn((n 1)(e o +2ir/3)+a n )+ srn((n-1 )(e o -2»/5)+« n )l} (X 3) 
It may be noted from (I 3) that two harmonic components 
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of phase current contribute a single harmonic component in 
the & axis current For example 5th and 7th harmonics m 
phase current contribute to 6th harmonic d axis current 
11th and 13th produces 12th harmonic d-axis current and so on 
The d axis current contain the harmonics of order 6 12 18 and 
so on Since the phase current harmonics are of 5 7 11 13 17 
etc The general expression for i^ n m terms of the Fourier 
Components of phase current as obtained from equation (13) is 

ha = / '57 T <( a n 1 + 3h + i) 2 + 0>n 1 + Vi )2) i (1 4) 

Similarly iq n can be derived as 

V " ^ {(a n-1 " a n+1 )2 + (lj n 1 Vl )2) + (I 5) 

where 

n « 6 12 18 
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appendix j 

EXPRESSION POR A<|> 

!, A 

| The non linear equation defining <{> (equation (5 ■&) of 

Chapter 5) can be written as 

J; 

0-j-j exp( BM) + C-j 2 exp( $ $ ) = 0 (J 1 ) 

where 

f 

n 

0^ » -/3/2 oos(<fr+6) (1+exp(-£ir/3))-sin(4’+6) exp($w/3) + 

8in(2ir/3 +<j>+$) + sin(ir/3+ ( l>+6) exp(-fto /3) 

C-j g ® /3/ cos(ot+6) (1+exp( 0 it/ 3)) sin(a+5)exp(- 3 2 n/ 3 ) 
sin( 2 ir /3 + ct + 6) sin( w /3 + a + 5) exp( £it/3) 

5 o 9 - p 0 = tan ^ (B^/Eg) = K 

Eg = V m + K ^ 2d P = tan 1 CX!,/^ ) 3 = V X 1 

On linearising equation (J 1) we get 

exp ( $a)( $ C-j-j Aa + A 0-j^ ) + exp ( B<|> ) ( 0 C^A^ + A 0^ 2 ) 35 0 (J2) 


and 




40 ii 

= c 22 (A<(>+Ae) 

(J 

3) 

A°1 2 

= C^j(Ao+A0) 

(J 

4) 

where 




° 22 

3/2 sin (<► + $) (1+exp( bV3))-cos(*+ 6) exp($ir/3) 

+ 



cos (2tr/3+ <fr +«)+oos(ir/3+^+«)exp( B"/3) 
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°33 = J 3/2 sin( a +<S )(1+exp( 0„/3))-oos( a+ S) exp ( e 2 „/j) . 
oos(2,/3+a + «) oos(n/3+a + 5) exp ( gj/j) 

Substituting (J 3) and (J 4) in (J 2) „ e get 

V* - c 5 a« + 0 6 U (J5) 

■where 

°4 ” °22 ex P(~ 6(* ♦)) 6 O-jg C 5 ■ 00^ exp( B(a <{>)) 

a 6 " °33 - °22 ex P< & ( a 4)) 

The equation defining 8 is given by 

E q tan e - £ d (J 6) 

On linearizing the above equation we get 


A6 . C ? A^ 2d + C 8 W 2q + 0 9 AY, 
•where 


m 


°7 - - 


b 4 k 




0, 


8 


E q g 

+ E n) 


B a 


C Q . — jH* Tjr 

7 . -nM 


J d T "V T Jj q'' ^ E d + V 

Eliminating A8 from equation (J 5) using (J 7) 

N 2d 


&<f> _ g 


T 


Ai p 


2q 


+ £3 Aa + £4 ^ ^ 


m 


where 


£ 


.T 


[ g-j “ 0 


°6 °7 


g 2 " 0 


C 6 °8 


°5 °6 °9 

g 3 " ^ g 4 * 


(J 7) 


(J 8) 



APPENDIX K 


SMALL SIGNAL MODEL OP AN INDUCTION MOTOR JITH 
SINUSOIDAL EXCITATION 

Let the source voltage v he defined by 

v « V„ sm(wt) 
ra 

so that 

V = 3 m \ 

Referring to Pig (5 2) of Chapter 5 current i is given by 

l 3 fEl '001 K ( 111213 + 3 Isal 

Z 1 

where 

2.j ** R,j + D 

Afber simplifio tion we get 

i a = m 1 K B, * 2q X 1 (V ffl + Kl' 24 )l/|zU 2 (B 

i . *^572 I K X 1 + K. » 2d )l 2 (! 

On linearizing (k 1) and (I 2) 
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where 


i z if 


x i E i 


l E i h 


5- he linear dynamical equation of the motor (ref eQuat on 
6 5)) is 


pfai|X + B^Y+C^u (K 4) 

where 

X T = lA¥ 2 a A* 2q AS] Y T (Ai a A lg] u- AT m 
and and 0^ are defined in equation (6 5) 

Substituting (K 3) in (K 4) we get 

p X = A X + C 1 u (K 5) 

where 

A = A^ + A 2 and A 2 = P 0 ) 
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